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We report and interpret HST /Space Telescope Imaging Spectrograph (STIS) 
long-slit observations of the optical and ultraviolet (UV) (1150 — 10270 A) 
emission-line spectra of the rapidly brightening Spot 1 on the equatorial ring 
of SN 1987A between 1997 September and 1999 October (days 3869 - 4606 after 
outburst). The emission is caused by radiative shocks created where the super- 
nova blast wave strikes dense gas protruding inward from the equatorial ring. We 
measure and tabulate line identifications, fluxes and, in some cases, line widths 
and shifts. We compute flux correction factors to account for substantial inter- 
stellar line absorption of several emission lines. 

Nebular analysis shows that optical emission lines come from a region of cool 
(T e ~ 10 4 K) and dense (n e ~ 10 6 cm -3 ) gas in the compressed photoionized 
layer behind the radiative shock. The observed line widths indicate that only 
shocks with shock velocities V s < 250 km s _1 have become radiative, while line 
ratios indicate that much of the emission must have come from yet slower (V s <^ 
135 km s" 1 ) shocks. Such slow shocks can be present only if the protrusion has 
atomic density n^,3x 10 4 cm" 3 , somewhat higher than that of the circumstellar 
ring. We are able to fit the UV fluxes with an idealized radiative shock model 
consisting of two shocks (V s = 135 and 250 km s" 1 ). The observed UV flux 
increase with time can be explained by the increase in shock surface areas as the 
blast wave overtakes more of the protrusion. The observed flux ratios of optical 
to highly-ionized UV lines are greater by a factor of ~ 2 — 3 than predictions 
from the radiative shock models and we discuss the possible causes. We also 
present models for the observed Ha line widths and profiles, which suggests that 
a chaotic flow exists in the photoionized regions of these shocks. We discuss what 
can be learned with future observations of all the spots present on the equatorial 
ring. 

Subject headings: supernovae: individual (SN 1987A) - supernova remnants - 
circumstellar matter 



1. Introduction 

Supernova 1987A (SN 1987A) in the Large Magellanic Cloud provides an unprecedented 
opportunity to observe the birth and development of a supernova remnant. International 
Ultraviolet Explorer (IUE) observations (Fransson et al. 1989; Sonneborn et al. 1997) found 
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narrow line emission about 80 days after the explosion 20 , demonstrating the presence of 
circumstellar gas around SN 1987A. Images taken with the Hubble Space Telescope (HST) 
showed that this gas consists of an equatorial inner ring (radius ~ 0.7 light year, n e ~ 
3 x 10 3 — 3 x 10 4 cm" 3 ) and two outer rings (~ 3 times the size of the inner ring, n e < 
2000 cm~ 3 ) tilted towards the observer at ~ 45° (Jakobsen et al. 1991; Wang 1991; Plait 
et al. 1995; Burrows et al. 1995; Lundqvist & Fransson 1996; Maran et al. 2000; Lundqvist 
& Sonneborn 2001). The circumstellar ring system was excited by the ionizing radiation from 
the supernova during the shock breakout (Lundqvist & Fransson 1991; Ensman & Burrows 
1992; Blinnikov et al. 2000). In the interacting winds model, the SN 1987A ring system is 
part of a bipolar shell around the supernova (Blondin & Lundqvist 1993; Martin & Arnett 
1995; Link et al. 2001). 

The first signal of ongoing interaction between the SN 1987A debris and the circumstellar 
gas was the rebirth of the supernova in X-ray (Beuermann, Brandt, & Pietsch 1994; Goren- 
stein, Hughes, & Tucker 1994; Hasinger, Aschenbach, & Triimper 1996) and radio (Staveley- 
Smith et al. 1992, 1993) wavelengths around day 1000. The size of the radio-emitting region 
indicated that the supernova debris expanded unimpeded at velocity j> 35, 000 km s _1 up 
to ~ day 1000 before slowing down to 3500 km s _1 by the interaction (Gaensler et al. 
1997; Manchester et al. 2001). The X-ray and radio observations have been explained as the 
interaction of the supernova ejecta with a rather dense (uh ~ 100 cm -3 ) H II region that 
separates the shocked stellar wind of the supernova progenitor from the denser gas of the 
inner ring (Chevalier & Dwarkadas 1995; Borkowski et al. 1997a). 

The "main event" of the birth of the Supernova Remnant 1987A (SNR 1987A) — the 
interaction between the supernova debris and the circumstellar rings — has been anticipated 
since the discovery of the circumstellar gas. Predictions of the time of the first contact ranged 
from 2003 (Luo & McCray 1991) to 1999 ± 3 (Luo, McCray, & Slavin 1994) to 2005 ± 3 
(Chevalier & Dwarkadas 1995). There have been previous model predictions of radiation 
from this impact in X-rays (Suzuki et al. 1993; Masai & Nomoto 1994; Borkowski et al. 
1997b) and in the UV/optical (Luo et al. 1994). The first definitive sign of impact between 
the supernova blast wave and the inner ring was detected in the 1997 April HST /Space 
Telescope Imaging Spectrograph (STIS) spectral images as a blueshifted (~ —250 km s _1 ) 
Ha feature at position angle (PA) th 30° of the ring (Sonneborn et al. 1998). Subsequent 
analysis of the HST /Wide Field and Planetary Camera 2 (WFPC2) images taken in 1997 
July showed that a point emission, located in projection at 88% of the distance to the ring, 
was increasing in brightness over a wide range of wavelengths (Garnavich et al. 1997, 2001) 



20 Thc time of core collapse of SN 1987A, 1987 February 23.316 (UT) (JD = 2,446,849.816), was determined 
by the 1MB and Kamiokande II neutrino detectors (Bionta et al. 1987; Hirata et al. 1987). 
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and could be traced back to as early as 1995 March (day 2932) (Lawrence et al. 2000a). The 
position of the brightening spot, located just inside the edge of the inner ring, suggested that 
this is the result of an interaction of the supernova blast wave with an inward protrusion 
of the ring. This brightening spot, named Spot 1, has increased in flux by a factor of ~ 10 
between 1996 and 1999 (Garnavich et al. 2001). With a number of new spots (~ 10 in 2000 
November) since then detected all around the inner circumstellar ring (Garnavich et al. 2000; 
Lawrence et al. 2000a), we are now witnessing the full birth of SNR 1987A. 

We present here HST /STIS UV and optical spectroscopy of Spot 1 taken by the Super- 
nova INtensive Study (SINS) collaboration up to day 4606 (1999 October 7). Results from 
an earlier (day 3869, 1997 September 27) STIS spectrum of Spot 1 have been presented in 
Michael et al. (2000). We describe the new observations and data reduction in §2 and report 
the results in §3. In §2.3 we describe a detailed method to determine the intrinsic fluxes 
and widths of a few UV lines, including the Si IV AA1394, 1403 and C IV AA1548, 1551 
doublets, which are strongly affected by interstellar line absorption along the line of sight to 
the supernova. 

Our working model (Figure 1) for Spot 1, as well as the other spots, is that it is caused by 
the impact of the supernova blast wave with a dense inward protrusion of the ring (Michael 
et al. 2000). When the blast wave overtakes such an obstacle, slower shocks are transmitted 
into it. Since a range of densities is present in the ring, we expect (§4.1) that the transmitted 
shocks will have a range of velocities (V s ~ 100 - 1000 km s" 1 ). Not all of these shocks are 
responsible for the observed UV/optical emission from Spot 1 though. While some UV and 
optical line emission is produced right at the shock front, a much larger amount is produced 
if the shocked gas undergoes thermal collapse, i.e. the shock becomes radiative (§4.2). The 
time it takes for a shock to make this transition increases with its speed and decreases with 
its pre-shock density. Once the post-shock gas collapses, a shock becomes an extremely 
efficient radiator of UV and optical emission lines (§4.3). Therefore, while the range of 
possible velocities present in the ring is large, we are only observing the shocks which are 
slow and/or dense enough to have become radiative. Michael et al. (2000) confirmed this 
general picture and found that the observed line widths and line intensity ratios indicated the 
emission was formed by radiative shocks in the velocity range ~ 100 — 250 km s -1 entering 
into a dense (n > 10 4 cm -3 ) gas. 

Nebular analysis of the observed emission lines of Spot 1 (§5.1) suggests that the emission 
comes from a region of high density (n e > 10 6 cm -3 ), This result confirms our picture that 
this emission comes from radiative shocks, which can compress the pre-shock gas by a factor 
;> 100. In §5.2, we model the UV line fluxes from Spot 1 between day 3869 and day 4587 
with the one-dimensional steady state shock code of Cox & Raymond (1985). The observed 
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UV fluxes are fitted well with models containing emission from two radiative shocks with 
V s = 135 and 250 km s^ 1 . We propose two scenarios to interpret the observed increase of UV 
flux with time of Spot 1. In the first scenario the density of the obstacle is low enough so that 
the cooling times of shocks entering the obstacle are comparable to the age of the spot. The 
increase in UV line fluxes are then due primarily to the aging of the shocks as they develop 
thermally collapsed layers. In the second scenario, the obstacle is dense enough so that all 
the shocks cool almost immediately. The increased observed fluxes are instead attributed to 
an increasing surface area of shock interaction. While the second scenario fits the observed 
fluxes better, we suspect that the actual light curves probably manifest a combination of both 
effects. In §5.3 we discuss the observed line profiles and present simulated Ha line profiles 
based on simple geometric models for the shock interaction. In §6 we discuss what we have 
learned by comparing results of simplified shock models with the spectral observations, and 
describe how future observations may elucidate some unsolved problems. We summarize the 
main results in §7. 



2. Observations and Data Reduction 

The STIS observations of Spot 1 in both optical and UV wavelengths obtained by the 
SINS team are listed in Table 1. SN 1987A is located in a densely populated region of the 
LMC and appears to belong to a loose, young cluster region (Panagia et al. 2000). Target 
acquisition was complicated by the stars present near the supernova, especially Star 3, a Be 
star of V ~ 16, at l'/63 away and PA = 118°, and Star 2, a B2 III star of V = 15.0, at 2'.'91 
away and PA = 318° from the supernova (Walborn et al. 1993; Scuderi et al. 1996). We 
decided to peak-up on the nearby star S2 (Walker & Suntzeff 1990) and offset the telescope 
to center the aperture on Spot 1. We measured the required offset from the WFPC2 images 
(Garnavich et al. 2001). Due to the uncertainties in this measurement, Spot 1 was located 
at 0'.'08 off the center of the slit for all observations taken before 1999 August. We reduced 
these data using the standard calibration files which assumed that the object was located at 
the center of the slit. We estimate that the offset from the center of the slit will cause us to 
underestimate the measured flux by ^ 5% for the 0'.'2 data, and by <^ 0.1% for the CK/5 data. 

With each grating setting, we took multiple (3 — 5) observations centered at dithered 
positions C'5 apart along the slit. Cosmic rays (in the case of optical data) and hot pixels 
(in optical and UV data) were removed simultaneously when the dithered raw images were 
combined with the CALSTIS software developed by the STIS Investigation Definition Team 
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at the Goddard Space Flight Center 21 . Previous narrow-slit STIS spectra processed by the 
SINS team with this software showed that the flux calibration of far-UV (G140L) and near- 
UV (G230L) data agree to < 2% for the overlapping region, while the agreement between 
near-UV (G230L) and optical (G430L) data is good to < 5% (Baron et al. 2000; Lentz et 
al. 2001). 

The location and orientation of the aperture positions are shown in Figure 2. For all 
but one of the observations, the slit was oriented (within ±5°) along the axis connecting 
the center of the SN 1987A debris and Spot 1, which is located at a PA= 29° on the inner 
ring (Garnavich et al. 2001). The only exception was the 1999 October G140L observation 
(Figure 2b), where the (52" x 0'/2) slit had a PA of 55° and did not pass through the 
center of the SN 1987A debris. In all observations, the Spot 1 spectrum overlapped with 
that from the adjacent segment of the inner circumstellar ring that was included within the 
STIS aperture. With an expansion velocity of 10.5 ± 0.3 km s" 1 (Cumming & Lundqvist 
1997; Crotts & Heathcote 2000), the ring was not resolved spectrally in any of the STIS 
observations reported here. In the optical data, the emission from the ring is the main 
source of background to the Spot 1 spectrum and will be discussed in §2.1 and §2.2. 

Garnavich et al. (2001) measured the width of Spot 1 in WFPC2 images up to 1999 
April 21 (day 4440) and showed that Spot 1 was unresolved in the data and was consistent 
with a point source at optical wavelengths. We compared the Full Width at Half-Maximum 
(FWHM) of Spot 1 emissions in our last STIS observations at 1999 September with those of 
point sources that were recorded in the data. We found that the FWHM of Spot 1 was 8 ±6% 
and 25 ± 9% larger than a point source in the far-UV and optical wavelengths, respectively. 
The latter result is consistent with measurements by Lawrence et al. (2000b) with the 2000 
May 1 (day 4816) STIS G750M spectroscopy, which suggests that Spot 1 is now moderately 
resolved in the HST data at optical wavelengths. 

2.1. Low Resolution Optical Observations 

We extracted the low resolution optical spectrum of Spot 1 from the STIS G430L and 
G750L two-dimensional spectral images. Portions of the 1999 September G430L and G750L 
data taken with the Cf!2 slit are shown in Figures 3(c) and 3(d), respectively. The horizontal 
streaks near the center in the figure are broad emission lines from the inner supernova 
debris, which have a FWHM velocity v fwhm « 2800 km s" 1 (Wang et al. 1996; Chugai et 
al. 1997). For these data, the lower section of the spectral image is the combined Spot 1 



CALSTIS Reference Guide, http://hires.gsfc.nasa.gov/stis/software/doc_manuals.html 
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and inner ring emission-line spectrum (hereinafter referred to as the Spot 1+ North- Ring, or 
Sl+NR, spectrum). At the kinematic resolution of the G430L and G750L grating settings 
(Ay « 300 — 550 km s^ 1 ), neither the emission lines from Spot 1 [vfwhm ~ 200 km s -1 , 
Michael et al. (2000)] or those from the ring (dfwhm — 10 km s _1 ) are resolved. The upper 
section of the spectral image is the emission-line spectrum of the segment of the inner ring 
subsection that is in the slit and directly opposite that of Spot 1 (hereinafter referred to as 
the South- Ring, or SR, spectrum). 

We measured the Sl+NR and SR spectra in the G430L and G750L grating settings by 
integrating the 7—9 rows of the image in which the emission-line spectra appeared. Emissions 
due to the diffuse LMC background in the Balmer lines, [O II] AA3727, 3729, and [O III] 
AA4959, 5007 are observed as images of the entire slit at those wavelengths in Figures 3(c) 
and 3(d). We subtracted the contribution of this diffuse emission from the Spot 1 spectrum 
by linear interpolation above and below the extracted rows. Since the subtracted LMC 
background was only a small fraction of the emission from Spot 1 (<2% for [O III] A5007), 
this subtraction did not significantly increase the uncertainty of the estimated line fluxes. 

We measured the flux of each emission line in the Sl+NR and SR spectra by fitting 
a Gaussian to the line profile, allowing the net flux, width, and central wavelength to vary 
independently for each line. We estimated the background level from a linear fit over a 45- 
pixel region of the spectrum centered on the line in question but excluding any emission-line 
features. We fitted the line flux by minimizing the total x 2 i n which the coefficients defining 
the line and the background were free parameters. Gaussian profiles gave satisfactory fits to 
all the line profiles. The reduced-^ 2 , or x 2 , (the total x 2 divided by the number of degrees 
of freedom), of our line fits were within the range 0.8 — 2.8, compared to the value of 1.0 
for a statistically good fit. We computed the flux of each line and its error from the best-fit 
parameters and their associated uncertainties. We adjusted the statistical errors of all the 
line fluxes so that x 2 — 1-0 f° r & U fits. In cases where two or more emission lines overlapped 
in wavelength, such as [O III] A5007 + He I A5016, and [Ar III] A7136 + [Fe II] AA7155, 7172, 
we fitted the emission features with multiple Gaussians with the additional constraint that 
their wavelength separations were the known differences of laboratory wavelengths. 

Only a subset (~ 1/3) of the emission lines observed in the Sl+NR spectrum also 
appeared in the SR spectrum. Therefore, we attributed entirely to Spot 1 the measured 
fluxes of emission lines that were seen in the Sl+NR spectrum but not in the SR one. 
To subtract the NR spectrum from the Sl+NR spectrum, we assumed that the fluxes of 
emission lines in the NR spectrum are equal to those in SR spectrum scaled by factors that 
are independent of time. This assumption is reasonable because the rate of flux decrease 
around the ring has been shown to be relatively constant around the ring (Plait et al. 1995; 
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Lundqvist & Sonneborn 2001). We estimated the scale factors by examining archival WFPC2 
emission-line images in Ha, [O III] A5007, and [N II] A6583 obtained in 1994 February and 
1994 September, before Spot 1 appeared. We measured flux ratios f(NR)/f(SR) of 1.2, 
1.1, and 1.3, respectively from the Ha, [O III] A5007, and [N II] A6583 images. We used 
the Ha scale factor for the Balmer lines. While we used the [O III] A5007 factor for the 
[O III] AA4363, 4959, 5007 lines, and likewise for the [N II] lines, we recognized the increase 
in systematic errors in the measured Spot 1 fluxes of [O III] A4363 and [N II] A5755 because 
these lines are more temperature sensitive than the other lines. For all remaining emission 
lines, such as [S II] and [Ne III], we assumed a scale factor of 1.2 ±0.1 in order to subtract 
the NR spectrum. 

In addition to the C'2 data sets described above, we obtained one observation of Spot 1 
with the 0'.'5 slit and the G750L grating setting. As before, we extracted the combined 
Sl+NR spectrum by integrating the 6 rows of the image where the emission lines appeared. 
For the emission lines produced predominantly by Spot 1, such as He I A6678 and [Ar III] 
A7135, we measured the SI fluxes by fitting the line emissions with single Gaussian profiles. 
For the lines where emission from both Sl+NR were apparent, such as [N II] A5755 and 
[O I] A6300, we estimated the contribution of the NR flux to the Sl+NR flux from a linear 
interpolation of the NR emission-line flux adjacent to Spot 1 in the slit. After subtracting 
this background, we fitted the remaining SI emission lines with Gaussian profiles. The 
uncertainties in the background ring flux estimated in this procedure resulted in larger 
systematic errors in the estimated Spot 1 line fluxes for this observation. 



2.2. Medium Resolution Optical Observations 

A medium resolution optical spectrum was taken on 1999 August 30 (UT, 4368.0 days 
since explosion) with a O'.'l slit and the G750M (6581) grating setting (6295-6867 A). With a 
spectral resolution of AV ~ 50 km s -1 , emission lines from Spot 1, with wfwhm ~ 200 km s" 1 
(Michael et al. 2000), were resolved in the data, while the emission from the unshocked inner 
ring, with fFWHM — 10 km s _1 , remained unresolved. Three observations of 7800 seconds 
each were taken at three parallel slit positions, pointed so that the middle slit position 
was centered on Spot 1 and the two other slit positions were immediately adjacent [cf. 
Figure 2(d)]. Therefore the observation covered a segment of the ring of length 0'.'3. 

With the crowded stellar field near SN 1987A, we did not execute the acquisition-peakup 
exposure for these O'.'l slit observations as suggested by the HST /STIS operation manual 
(Leitherer et al. 2000). Instead all three adjacent slit positions were placed on Spot 1 by 
blind offsets. Therefore we cannot apply the standard pipeline data reduction procedures to 
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process these data. To determine the fluxes of emission lines from Spot 1, we first removed 
the wavelength-dependent aperture throughput correction function in the pipeline data for 
each O'.'l slit spectrum. We then summed the fluxes from the three OVl slit positions, and 
multiplied the total flux by a new aperture correction function for an equivalent 0'.'3 slit, 
calculated by interpolating the pipeline corrections for the O'.'l, 0'.'2, and (X'5 slits. For Ha, 
the aperture correction led to a ~ 20% decrease of flux over the sum of the fluxes measured 
in the three slit positions. We calculated the corresponding la errors by combining the 
individual errors in quadrature. 

A section of the spectral image from the middle slit position is shown in Figure 3(e). 
Again, the central horizontal streak is emission from the SN 1987A debris. Emissions from 
the inner ring at the two positions where the ring intersected with the slit aperture were 
observed in [O I] AA6300, 6364, Ha, [N II] AA6548, 6583, and [S II] AA6717, 6731. In 
the lower section of the spectral image, emission from the stationary inner ring overlapped 
with the broadened emission from Spot 1. Again, the major source of contamination the 
Spot 1 spectrum is the emission from the inner circumstellar ring within the O'.'l slit. As 
described by Michael et al. (2000), we fit all the Spot 1 emission features with Gaussian 
functions. Emission from the stationary ring dominated the spectral profile near zero velocity 
(~ ±50 km s" 1 ) and was excluded from the fit. The majority of the line profiles could be 
fitted well with Gaussian profiles, such as the fit to the [N II] A6583 line emission shown in 
Figure 4. The sole exception was the Ha line profile, where the signal was strong enough to 
show noticeable departures from a Gaussian profile, as we shall discuss further in §5.3. 



2.3. Low Resolution UV Observations 

We obtained low resolution UV spectra of Spot 1 with the G140L and G230L grating 
settings. Michael et al. (2000) have already presented results from the first G140L far-UV 
observations in 1997 September 27 taken with the 0'.'5 slit. We detected emission lines from 
Spot 1 in N V AA1239, 1243, Si IV AA1394, 1403, O IV] A1400, N IV] AA1483, 1487, C IV 
AA1548, 1551, and He II A1640. We detected the same set of emission lines in 1999 February 
27, observing again with the 0'.'5 slit. In 1999 October 7, observing with the 0'.'2 slit, we also 
detected the C II A1335 multiplet, [Ne IV] A1602, and O III] AA1661, 1666. Figure 3(a) shows 
a section of the spectral image from this observation. Radiation from Spot 1 is evident in the 
lower portion of the image whereas only faint line emission from the inner circumstellar ring 
can be seen in the upper half of the displayed image. The broad (~ ±15,000 km s _1 ) Lya 
radiation comes from the reverse shock from the interaction between the supernova debris 
and the H II region located inside of the equatorial ring (Sonneborn et al. 1998; Michael 
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et al. 1998a,b). Fluxes of UV line emission from the inner ring are much less than those 
from Spot 1 and make a negligible contribution to the measured fluxes. This is also the case 
for the near-UV emission lines measured in the 1999 September 17 G230L 0'/2 observation, 
shown in Figure 3(b). 

We measured the far-UV and near-UV spectra of Spot 1 from the G140L and G230L 
data, respectively, by procedures similar to those we described in §2.1. We fitted emission 
lines with Gaussian profiles and linear backgrounds, except for the broad Lya emission 
underlying the N V AA1239, 1243 doublet which we fitted with a quadratic function. We 
fitted the two components of close doublets such as N V AA1239, 1243, C IV AA1548, 1551, 
O III] AA1661, 1666, N II] AA2139, 2143, and Mg II AA2796, 2803 with Gaussians constrained 
to have fixed doublet separations, identical widths, and line ratios dictated by their oscillator 
strengths. At the resolution of the G140L grating setting, the Si IV A1403 emission of the 
Si IV AA1394, 1403 doublet is blended with the O IV] A1400 multiplets. We fit the combined 
Si IV and O IV] feature with multiple Gaussians, requiring the Si IV doublet to meet the 
same constraints as the other close doublets. 

The observed fluxes of a few UV lines, such as Si IV AA1394, 1403 and C IV AA1548, 
1551, were reduced by interstellar line absorption. We describe our procedure for correcting 
for this absorption in §3.3 below. 

3. Data 

3.1. Optical Emission Line Fluxes 

Table 2 lists the measured Spot 1 optical emission line fluxes, including previously 
published results from the 1998 March 7 (day 4030) data by Michael et al. (2000) and la 
upper flux limits for the [O II] AA3726, 3729, and [N I] AA5198, 5200 doublets. The la errors 
tabulated are only statistical errors. Systematic effects, such as fringing for the G750L data 
towards the near-IR region (A > 8500 A, Leitherer et al. 2000), might contribute additional 
uncertainties to the measured fluxes. 

The tabulated fluxes have been dereddened with E(B — V) of 0.16 (Fitzpatrick & Wal- 
born 1990) and the extinction law of Cardelli et al. (1989) with an assumed Ry of 3.1. In the 
optical band the differences between the LMC extinction law and the Galactic law are neg- 
ligible at low color excess (Fitzpatrick 1999). The interstellar extinction correction applied 
is listed in the last column of Table 2. 

Spot 1 was observed in many neutral and lowly ionized species in the optical wavelengths. 
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We did not detect any coronal lines such as [Fe X] A6375. Most of the line fluxes increased 
with time at a rate comparable to that measured from WFPC2 photometry (Garnavich et 
al. 2001). At the low spectral resolution of both the G430L and G750L observing modes 
(R = A/AA = 530 — 1040), definitive line identification remained a problem, especially for 
the [Fe II] emission lines. The Fe line identifications in Table 2 are based upon the modeling 
of the Spot 1 Fe lines in Pun et al. (2002, in preparation). Moreover, several lines were 
blended. Table 2 lists possible contributing species and, in the case of [Fe II] lines, different 
multiplets. In contrast, line blending is not a problem in the medium resolution G750M 
{R « 6000) data. 

Fluxes in [O I] AA6300, 6363 and He I A6678 were measured with the low resolution 
G750L and medium resolution G750M gratings in 1999 September within 17 days of each 
other. The two measurements agreed within uncertainties for the [O I] doublet. The two 
results differed by ~ 50% for the low S/N He I A6678 data, but were also within the noise 
level. This difference is probably a good indication of the detection limits of such faint lines. 

3.2. Optical Emission Line Widths 

For the medium resolution G750M observations, apart from the line fluxes, we were also 
able to measure the peak emission velocities (Vo) and the widths (Vfwhm) of the emission 
lines from the profile fits. Table 3 lists the results. The peak emission velocity measurements 
have been adjusted for the SN 1987A heliocentric velocity of +286 km s _1 (Wampler et al. 
1989). The 1998 March / .'2 G750M results have also been adjusted for the off-center position 
of Spot 1 within the aperture (cf. §2). The main uncertainties in the measurements of both Vq 
and Vfwhm are due to the contributions by emission from the stationary circumstellar ring, 
which dominated the emission by Spot 1 near zero velocity for many species (cf. Figure 4). 
The errors due to this contribution are generally smaller for the 1999 August O'.'l observations 
than the 1998 March 0'.'2 ones. 

For all emission lines, the line centroids from Spot 1 were blueshifted, with peak velocity 
Vq lying within the range —40 to —10 km s -1 . This result is consistent with the overall 
physical picture in which Spot 1 is located on the near side of the equatorial ring (Sonneborn 
et al. 1998) and the shock entering Spot 1 is moving towards the observer. We found no 
evidence that the peak velocity V of Spot 1 varied with time. 

Most emission lines had widths (FWHM) within the range ~ 150 — 180 km s -1 , except 
[N II] A6583 and Ha, which had a slightly greater width, Vfwhm ^ 200 km s -1 . We detected 
no measurable change with time of the line widths except for Ha and [N II] A6548. For Ha, 
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the emission profile from the second observation in 1999 August could no longer be fit well by 
a single Gaussian (to be discussed further in §5.3). The width of the [N II] A6548 emission, 
measured in 1999 August was 1.5 times greater than that measured in 1998. However, we 
are inclined to attribute this increase to systematic error, since we detected no such increase 
in the other [N II] component at 6583 A, where the line fluxes were ~ 3 times stronger. 

3.3. UV Emission Lines 

3.3.1. Interstellar Line Absorption 

Near the time of outburst, interstellar absorption lines of C II A1335 multiplet, Si IV 
AA1394, 1403, C IV AA1548, 1551, and Mg II AA2796, 2803 were observed in the UV con- 
tinuum of SN 1987A with IUE operating in the high resolution (AV ~ 30 km s^ 1 ) echelle 
mode (Blades et al. 1988; Welty et al. 1999). For each line, the dominant absorption compo- 
nent was centered at +281 km s -1 , near the SN 1987A heliocentric velocity of +286 km s _1 
(Wampler et al. 1989), and had a FWHM of ~ 80 km s -1 . Therefore, narrow emission lines 
(vfwhm — 10 km s _1 ) from the inner circumstellar ring from these species are completely 
blocked by the interstellar absorption (Fransson et al. 1989), as demonstrated by the absence 
of Si IV and C IV emission lines from the ring in Figure 3(a). 

Emission lines from Spot 1, with Vfwhm ~ 200 km s _1 , are not totally blocked by these 
interstellar absorption lines, as correctly predicted by Luo et al. (1994). However, the line 
profiles are altered and the observed fluxes are reduced substantially. Therefore, we must 
correct the measured C II, Si IV, C IV, and Mg II emission line fluxes to account for this 
absorption. The appropriate correction factors depend on the profile shapes of both the 
Spot 1 emission and the intervening absorption. To make this correction, we assumed that 
the UV emission lines from Spot 1 had Gaussian profiles with the same parameters as the 
optical lines as measured in the medium resolution optical data (§2.2). We then estimated 
the amount of flux reduction in the UV emission lines by multiplying the assumed Gaussian 
profiles by the absorption profiles seen in the IUE data. The corresponding flux correction 
factors for the various emission lines are shown in Figure 5(a). The correction factor is 
greatest for the Mg II AA2796, 2803 doublet, where the interstellar absorption is almost 
completely saturated between —50 and +300 km s -1 . 

The flux correction factor due to interstellar absorption is sensitive to the assumed 
widths (Vfwhm) an d peak (Vo) velocities of assumed line profiles from Spot 1. Figure 5(b) 
illustrates this dependence for the important case of C IV AA1548, 1551. We see that the 
flux correction factor increases moderately with decreasing Vfwhm for Vfwhm > 200 km s _1 , 
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but becomes very sensitive to both Vfwhm and Vq for Vfwhm < 150 km s 1 . 

As we shall discuss below in §5.3, the peak velocities and widths of the emission lines 
from Spot 1 depend on the detailed geometry and hydrodynamics of the shocks entering the 
spot, which are unknown. It is not obvious that the optical and UV emission lines should 
have the same peak velocities and widths. However, in the plane-parallel shock model that 
we describe in §4, we found that the peak velocities of emission from Spot 1 were almost 
identical in both the UV and optical wavelengths. Therefore, we used the measured peak 
Spot 1 velocities from the medium resolution optical observations, Vq = —30 ± 15 km s _1 
(cf. Table 3), to estimate the flux correction factors to account for interstellar absorption of 
the UV emission lines. 



3.3.2. Widths of the UV Emission Lines 

The far-UV emission lines from Spot 1 are poorly resolved at the spectral resolution of 
the G140L grating (AV ~ 250 km s _1 at 1500 A). For these lines, we attempted to establish 
the relation between the actual and observed widths through Monte Carlo simulations. For 
the simulations, we assumed that the intrinsic Spot 1 emission lines had (i) a Gaussian shape 
with FWHM as a variable parameter, (ii) emission peaked at —30 km s _1 , and (iii) flux ratio 
of the doublets dictated by their oscillator strengths. 

The dotted curve in Figure 6(a) shows a model C IV A A 1548, 1551 profile, assuming 
Vfwhm = 150 km s _1 and line flux ratio J(1548)//(1551) = 2:1. The solid curve shows 
a subsection of the high resolution IUE spectrum with the absorption line profile for the 
C IV doublet. For each input FWHM, we multiplied the model input emission line by the 
measured IUE absorption spectrum. Figure 6(b) shows a typical result. We convolved this 
profile with the line-spread function (LSF) of the detector for the slit aperture of the data 
set. The LSF for each emission line was interpolated from the measured LSFs at 1200 A, 
1500 A, and 1600 A. The LSF typically has a narrow peak (~ 1.5 pixel) and a broad wing 
that extends to ~ 10 pixels on either side of the peak (Leitherer et al. 2000). Figure 6(c) 
shows a typical example of the LSF for the G140L grating at 1550 A with the 0'/2 aperture. 

For each emission line, we normalized the resulting convolved profile to the photon 
counts for each observation. The thick solid curve in Figure 6(d) shows the normalized 
model C IV line profile for the 1999 October observation. We constructed a simulated 
observed profile (the thin solid curve in Figure 6d) by applying random noise to the model 
profile and sampling it at the spectral resolution of the G140L grating. Finally, we fitted 
the simulated profile with Gaussians in the same way as the real data. The dotted curve in 
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Figure 6(d) shows such a fit. 

For each input model FWHM velocity, we ran 10000 simulations and generated an array 
of the corresponding observed line widths. Figure 7 shows the median and the 68.3% (la) 
upper and lower limits of the array plotted against the model input widths of the various 
observed emission lines. We use these results to convert the observed line widths from the 
STIS data, shown as horizontal lines in Figure 7, to the intrinsic widths of the lines emitted 
by Spot 1. 

We did not attempt to model the near-UV G230L data in this way because the spectral 
resolution for this grating setting (AV 450 km s _1 at 2500 A) was too low for such a 
procedure to yield useful results. 

For all emission lines except the N V doublet, we found a la upper limit of FWHM < 
300 km s^ 1 , consistent with the measurements from the medium resolution optical data. 
The lower limit to the assumed FWHM is important because the flux correction factors for 
all UV emission lines are sensitive to the input line widths for FWHM < 150 km s _1 . We 
took this lower limit to be the same as the measured lower limit for the optical lines, that 
is, FWHM > 100 km s -1 . 

3.3.3. UV Line Fluxes and Line Widths 

Table 4 lists the fluxes of UV lines from Spot 1 inferred from the Gaussian line fits (cor- 
rected for extinction) and the intrinsic line widths derived from the Monte Carlo simulations. 
The tabulated fluxes for the C II, Si IV, C IV, and Mg II doublet have been corrected for 
interstellar line absorption assuming an intrinsic Spot 1 line width Vfwhm = 150^5q° km s _1 
with peak emission at Vq = —30 ± 15 km s _1 . The flux correction factors applied for these 
lines are also listed in Table 4. We also corrected the fluxes for interstellar extinction by 
assuming E(B - V) LMC = 0.06, E(B - V) Galactic = 0.10, and R v = 3.1. At UV wavelengths, 
the choice of extinction function is important because the correction is substantial and is 
known to vary from place to place (cf. Pun et al. 1995). We used the 30 Doradus extinction 
function of Fitzpatrick (1986) for the LMC component, and the Seaton (1979) function for 
the Galactic component. 

The far-UV G140L data had slightly higher kinematic resolution (AV « 250 km s _1 ) 
than the low resolution optical G430L and G750L data. However, the individual components 
of C II A1335 and O IV]A1400 multiplets remained unresolved in the data. There were more 
uncertainties with line identifications in the lower resolution (AV ~ 300 — 650 km s _1 ) 
near-UV G230L data, such as the unidentified emission features near 2737 A and 2746 A. 
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On the other hand, we identified emission features near 2324 A and 2334 A as the C II 
A2325 multiplet and Si II A2335, respectively. We ruled out the alternative identification 
of [O III] AA2322, 2332 because the observed line ratios, 7(2324 A)/7(2334 A) « 5 and 
7(2324 A)/7(4363 A) « 7, were much different from the theoretical [O III] line ratios of 
~ 280 and 0.12, respectively. These [O III] line ratios are determined only by the atomic 
transition probabilities and are independent of the excitation model. 

The fluxes of all far-UV lines from Spot 1 increased with time during the three obser- 
vations taken from 1997 September (day 3869) to 1999 October (day 4606). The rate of 
increase differed for features from different ions, ranging from 7(3869 d)/7(4596 d) = 5.3 
for N V A1240 to 1.9 for C IV A1550. We will discuss the time dependence of the UV 
emission-line fluxes in §5.2. 

4. Interpretation 

4.1. Impact Hydrodynamics 

In our working model for Spot 1, we assume that the UV and optical emission lines 
observed are caused by radiative shocks that develop where the supernova blast wave strikes 
dense gas protruding inward from the circumstellar ring. As we shall show, the spectrum 
and profiles of the emission lines from such shocks are sensitive to the density distribution 
and geometry of this protrusion, which are probably quite complex and cannot be resolved 
even with the HST. Our approach here therefore is to illustrate the salient physics of the 
spectrum formation with a few idealized "toy models," which we believe will guide us toward 
a better understanding of the properties of a more realistic model. 

Following the previous work of Luo et al. (1994) and Borkowski et al. (1997b), we show 
in Figure 8 hydrodynamic simulations for two models of a fast shock overtaking a dense gas 
cloud. In each model, we assume that a fast plane-parallel blast wave traveling through a 
uniform medium of relatively low density overtakes a cloud of substantially greater uniform 
density. The cloud boundary is approximated as a density discontinuity. In each simulation, 
the blast wave drives a transmitted shock into the cloud, while a reflected shock travels 
backwards towards the interior of the remnant. As the blast wave overtakes the obstacle, 
the surface area of the transmitted shock increases. The transmitted shock propagates with 
a range of velocities depending on shape and density of the obstacle. 

The simulations are calculated using the hydrodynamic code VH-1 (Strickland & Blondin 
1995), which is based on the piecewise parabolic method of Colella & Woodward (1984). Ra- 
diative cooling is included in the code using an operator splitting technique. We have used 
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a non-equilibrium ionization cooling curve calculated with the plane-parallel shock code de- 
scribed in §4.2 with abundances typical of the equatorial ring. At the resolution of the 
simulations, we are not always able to resolve the cooling time scale of the shocks. However, 
this computational limitation does not seriously affect the behavior of the overall hydrody- 
namics of the interaction. Moreover, we have modeled the blast wave as a single planar shock 
rather than using the full double-shock structure present in the remnant. We will comment 
on the effects of ignoring the full double-shock structure in §6. 

The two scenarios depicted in Figure 8 show two different behaviors for the development 
of radiative shocks in the obstacle. The important parameter distinguishing their behavior 
is t cool /t cross , where t coo i is the typical cooling time for the transmitted shocks, and t cross is 
the time it takes for the fast shock to cross the obstacle. As we will discuss below in §5.2.2, 
the interpretation of the observed light curves of lines emitted from the shocks depends on 
which behavior is dominant in Spot 1. 

The simulation on the left of Figure 8 (Scenario 1) is one in which t coo i is comparable 
to tcross-, so that not all of the shocks have become radiative. For this simulation we assume 
that the obstacle is a spherical cloud of density p — 10 4 amu cm -3 and diameter d spo t = 
4x 10 16 cm. This geometry allows the the shock to be driven into the back side of the obstacle. 
Even after the blast wave has completely overtaken the obstacle, the shocks transmitted into 
the front-end of the obstacle have not yet undergone thermal collapse and are non-radiative 
(NR) because of the high impact velocity. On the other hand, the shocks transmitted into the 
sides and back of the obstacle have lower velocities owing to the oblique incidence of the blast 
wave. The shocked gas behind these parts of the transmitted shock has a lower temperature 
and a shorter radiative cooling time. A dense layer of cooled gas (R) is developed as a 
consequence (to be discussed further in §5.2). 

Scenario 2 in Figure 8 represents an example in which the cooling time for the trans- 
mitted shock is much shorter than the time scale for the blast wave to cross the cloud 
tcooi "C t cross , regardless of incidence angle of the blast wave. We assume that the obstacle is 
an elongated protrusion with a higher density, p — 10 5 amu cm -3 . With the high density 
of the obstacle, all the transmitted shocks undergo thermal collapse soon after impact. 

We recognize that these idealized models in Figure 8 do not represent the true shape 
and density distribution of the obstacle. However, given the limited observations at hand, we 
believe it is more fruitful to explore how well we can fit the data with a few idealized models or 
combinations thereof, rather than to explore fits of more complicated hydrodynamic models 
(see §6). 
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4-1.1. Transmitted Shock Velocities 

The driving pressure Pb behind a hypersonic blast wave propagating with velocity V& 
through a medium of density p is given by Pb = 3/4pV b 2 . Immediately after the blast wave 
passes the surface of an obstacle, a transmitted shock propagates into the obstacle with 
velocity 

V s = (Ps/p ) 1/2 , 

where P s is the pressure at the surface of the obstacle and p is the pre-shock density of the 
obstacle. The value of P s at a given location on the cloud surface rises shapely as the blast 
wave first strikes it, then decreases as the reflected shock propagates away from it. Over a 
time scale comparable to that for the blast wave to overtake the obstacle, P s decreases by a 
factor of ~ 1.5 — 2 (Borkowski et al. 1997b). Immediately after the passage of the blast wave, 
the ratio P s /Pb at any given point on the cloud surface depends on two parameters, first, the 
obliquity angle 9 between the direction of the blast and the inward normal to the surface of 
obstacle, and second, the pre-shock density ratio 5 = Po/phii between the obstacle and the 
ambient H II region. Figure 9 illustrates the dependence on obliquity for a blast wave with 
a density contrast 5 70 impacting on a spherical cloud. At the point of first contact, i.e., 
9 = 0, the driving pressure is at its maximum value and P s /Pb ~ 4. As the obliquity of the 
impact increases, the driving pressure decreases. With a factor of ~ 10 decrease in pressure 
due to obliquity, the initial velocity of the transmitted shock on the edges of the protrusion 
is only ~ 30% of that at the tip of the protrusion. 

We estimate the pressure behind the blast wave, Pb, from the relation P b = 3/ApV b 2 . 
By fitting the radio remnant of SN 1987A from year 1992 to 1995 with a shell model, 
Gaensler et al. (1997) measured the velocity of the blast wave to be Vj, = 2800 ±400 km s _1 . 
Using more recent radio observations, Manchester et al. (2001) updated the value of Vb to 
be 3500 ± 100 km s _1 . However, they also noticed that in data taken since ~ 1998, the 
radio remnant can no longer be fit well with the simple shell model but instead with a 
combined shell and multiple point-sources model. Manchester et al. (2001) estimated that 
this leads to a 30% uncertainty in the measured expansion velocity. On the other hand, by 
modeling the observed X-ray emission from the blast wave before any spots appeared on the 
circumstellar ring, Borkowski et al. (1997b) found a slightly higher velocity for the blast wave, 
< Vb >= 4100 km s _1 , and a density of the H II region inside the ring, p H n = 150 amu cm~ 3 . 
These results are consistent with the estimates of Chevalier & Dwarkadas (1995) and the 
upper limit pnn < 240 amu cm~ 3 determined by Lundqvist (1999). We estimate the blast- 
wave pressure to be within the range Pb — (1 - 5) x 10~ 5 dyne cm -2 . 

The pre-shock density, p , of the obstacle is even more uncertain. The density of the ring 
has been determined from the rate of fading of the optical and UV emission lines (Lundqvist 
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& Fransson 1996; Sonneborn et al. 1997). At the time of emission maximum (~ day 350), the 
radiation was dominated by gas of relatively high density, po ~ 5 x 10 4 amu cm" 3 . However, 
the emission from the higher density gas faded rapidly and the emission at later times was 
dominated by gas of lower density, po ~ 10 4 amu cm" 3 (Maran et al. 2000). The gas in 
Spot 1 might have even higher density than the value derived from optical and UV lines. 
The fact that Spot 1 lies inside the inner circumstellar ring might be due to the fact that it 
resisted ablation from the ionizing radiation and stellar wind of the progenitor as a result of 
its enhanced density. 

We summarize the possible range of transmitted shock velocities in Figure 10. The 
shaded gray region shows the range of transmitted shock velocity, V s , as a function of the 
pre-shock density of the obstacle, p , for our best estimate of the blast-wave pressure (V b = 
3500 km s _1 , p Hn = 150 amu cm -3 ). The upper boundary is the velocity of the transmitted 
shock at the tip of the protrusion and the lower boundary is the velocity at the side. We 
also show in Figure 10 the corresponding shock velocity range for our high estimate of 
the blast-wave pressure (Vb = 4100 km s _1 ,phii = 250 amu cm" 3 , dashed), and our low 
pressure estimate (V& = 2800 km s _1 ,phii = 100 amu cm" 3 , dotted). Assuming a pre-shock 
obstacle density of (1 — 5) x 10 4 amu cm" 3 , it is apparent that a range of shock velocities 
(~ 100 — 1000 km s" 1 ) can be present in the obstacle. 



4.2. Shock Structure 

In this section we review the relevant physics of shock fronts that we need to interpret 
the emission-line spectrum of Spot 1. We use the 1991 version of the Raymond shock code 
(Cox & Raymond 1985) for illustration. The code calculates the non-equilibrium ionization 
and excitation of the post-shock flow in a one-dimensional steady state shock for any given 
input parameters such as shock velocity and pre-shock densities. For any specified electron 
(T e ) and ion (Tj) temperatures at the shock front, the code calculates their equilibration 
through Coulomb collisions. We assume a ratio of T e /Tj = 0.2 for our models. The code 
calculates the upstream equilibrium preionization and the downstream photoionization of 
the gas and derives the post-shock density, temperature, and ionization structure in the gas. 
The code also calculates the local emissivity throughout the shock and the integrated fluxes 
over the entire column behind the shock front, including continuum and line emission from 
both allowed and forbidden transitions of H, He, C, N, O, Ne, Mg, Si, S, Ar, Ca, Fe, and Ni. 

We adopted the LMC abundances measured by Russell & Dopita (1992) in all elements 
except for He, C, N, and O, where the ring abundances derived by Lundqvist & Fransson 
(1996) are used instead. The other exception is Si where the abundance of Welty et al. 
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(1999) is used because of the large uncertainty in the measurements by Russell & Dopita 
(1992). We call this set of values the "Ring" abundance (H : He : C : N : : Ne : Mg : Si : 
S : Ar : Ca : Fe : Ni = 1 : 0.25 : 3.24 xl(T 5 : 1.82 xl(T 4 : 1.58 xlO" 4 : 4.07 xlO" 5 : 2.95 
xlO- 5 : 2.51 xl0~ 5 : 5.01 xl0~ 6 : 1.95xl0- 6 : 7.76 xlO" 7 : 1.70 xl0~ 5 : 1.10 xlO" 6 ). 

We show the post-shock temperature and density structure for protons and electrons 
in our model shock in the upper panel of Figure 11. We assume a shock velocity of V s = 
250 km s" 1 , characteristic of the shocks producing the radiation seen from Spot 1 (Michael et 
al. 2000). The shock has a high Mach number (M = V s /c s 3> 1), and compresses the obstacle 
gas by a factor of ~ 4 at the shock front, providing that the magnetic field is negligible. 
The temperature of the ions after crossing the shock front is given as 7$ = 3mjV^ 2 /16k, 
where rrii is the mass of the ion and k is the Boltzmann's constant. This implies that each 
ionic species will have a different post-shock temperature, provided that there are processes 
to completely thermalize the post-shock distribution functions. For collisional shocks, the 
post-shock electron temperature is orders of magnitude lower than the ion temperatures 
(Zeldovich & Raizer 1967). On the other hand, in collisionless shocks, plasma turbulence 
(e.g., Cargill & Papadopoulos 1988) can partially equilibrate the post-shock electron and 
ion temperatures. Eventually, the temperatures will be equilibrated by Coulomb collisions 
to a mean shock temperature T s = 3fim p Vg /16k, where /i — 4/(8 — 5Y) is the mean atomic 
weight per particle and Y is the mass fraction of helium. 

In the region (the "ionization zone") immediately behind the shock front with a charac- 
teristic length of a few ionization lengths, the atoms are collisionally ionized. Since radiative 
processes in this region remove a negligible fraction of the thermal energy of the hot plasma, 
this part of the shock is called the "adiabatic" or "non-radiative" zone. Given sufficient time 
the radiative losses will remove the thermal energy and cause a runaway thermal collapse of 
the shocked gas. In order to maintain pressure equilibrium across the shock in this "cooling 
region," an increase in density will accompany with a temperature decrease in the region. 
Radiative cooling in the gas continues until its temperature (T pa 10 4 K) is too low for 
the excitation of UV resonance lines which are responsible for the rapid cooling. Moreover, 
once the gas starts to recombine, it quickly becomes optically thick to the ionizing radiation 
(mostly extreme UV lines) produced upstream in the cooling zone. Roughly half of this 
ionizing radiation propagates downstream and is reprocessed into optical emission lines in 
this "photoionization zone," where the gas is maintained at a temperature T c pa 10 4 K as a 
result of a balance between photoionization heating and radiative cooling. 

The density of the gas in the photoionization zone is given by p c = 16p T s /3T c . For 
a V s = 250 km s _1 shock, where T s pa 10 6 K, the gas in the photoionization zone can be 
compressed by a factor pa 550. However, a significant magnetic field entrained in the gas 
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may mitigate this compression (Raymond 1979). 

Shocks that have developed cooling and photoionization zones are called "radiative 
shocks." A shock will be radiative if it has propagated for a characteristic cooling time, 
t coo i, which is defined as the time required for the shocked gas to cool from is post-shock 
temperature to 10 4 K. Otherwise, we call it a "non-radiative shock." We calculate the 
cooling time T coo \ for shocks in our models for input shock velocities V s = 100 — 600 km s _1 
and plot the results in Figure 12. We fit a power law to these results and obtain the relation 



where p is the pre-shock density. Therefore, for the range of shock velocities and pre- 
shock densities expected in Spot 1, which was less than 4 years old when our observations 
were made, some shocks may have already become radiative while others may remain non- 
radiative. 



After entering the shock front, atoms are collisionally ionized until they come into equi- 
librium with the post-shock gas. These ions are collisionally excited and emit line radiation, 
some of which from lower ionization stages than the final ionization state of a species. This 
ionization process in the ionization zone is not in equilibrium. Emission from this region 
lead to the detection of "Balmer filaments" in other supernova remnants (Hester et al. 1986; 
Long & Blair 1990). The Balmer emission comes from the shock front, where the H atoms 
are ionized. 

When ions reach equilibrium with the post-shock gas, the dominant line emission are 
in far UV or soft X-ray wavelengths, depending on the post-shock temperature. Most of the 
UV and optical line emissions are produced after the shock starts its thermal collapse. While 
the observed high-ionization UV lines are produced in the cooling region, the low-ionization 
UV and optical lines are produced in the photoionization zone, which is a dense H II region 
that is illuminated by the harder ionizing spectrum created upstream. 

We plot the integrated surface emissivities of N V A1240, C IV A1550, [N II] AA6548, 
6584, and Ha of our shock model described above in §4.2 as functions of downstream column 
density in the lower panel of Figure 11. Most of these lines (> 95%) are emitted during or 
after the thermal collapse of the shock, rather than in the ionization zone directly behind the 
shock front. Therefore, the resulting spectrum and magnitude of the emission from a shock 
depend greatly on whether the shock is radiative or non-radiative and the cooling time of a 




4.3. 



Shock Emission 
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shock is a good indication of the time it takes for a shock to "light up" . 
The luminosity of a line produced by a shock in ergs s _1 is given by 

L = ^V(y s , Po,t s hock)poV s 3 A s , 

where po is the pre-shock density, V s is the shock speed, t s h oc k is the age (time since first 
encounter) of the shock, and A s is the surface area that the shock covers. The function rj 
represents the fractional efficiency for a shock to convert its thermal energy into emission 
for a given line. As shown in Figure 11, 77 can be represented as a step function which 
turns on at t s h oc k = t CO oi- We show the relation between the emission efficiency 77 and shock 
velocity V s for several emission lines in Figure 13. The thick lines in Figure 13 represent 
shocks that have completely cooled, that is, t s h oc k 3> t coo i, while the thin lines are shocks 
with r = 3 x (2 x 10 4 amu cm~ 3 / p ) yrs. Shocks with velocities V s > 250 km s _1 have not 
yet cooled by t s h oc k = 3 yrs and therefore their structures are truncated when compared to 
the fully developed radiative shocks. This shows again that radiative shocks are far more 
efficient radiators than non-radiative shocks. Therefore, while non-radiative shocks may be 
present in the protrusion, their net contribution to the observed UV and optical emission 
from Spot 1 will be negligible. 

For emission lines such as N V A 1240 (formed in the cooling region), and Ha (formed 
primarily in the photoionization zone), the line efficiency 77 has little dependence on the 
pre-shock density p . In contrast, forbidden lines are subject to collisional supression and 
therefore their emissivities are sensitive to po- The effect of supression is illustrated in Fig- 
ure 13 by the [N II] AA6548, 6584 lines, which have a critical electron density of ~ 10 5 cm -3 
(Osterbrock 1989). Effects from collisional supression increases as the shock velocity increases 
because faster radiative shocks generate more compression. Figure 13 also shows that line 
emissions such as N V A1240 can be produced only if the temperature of the shocked gas is 
high enough for that ion to be produced, as previously discussed in Michael et al. (2000). 
Above that threshold, the total line emissivity increases linearly with the shock velocity. 
Permitted emission lines formed in the photoionization zone have a stronger dependence on 
shock speed, since they are generated by the reprocessing of the ionizing photons produced 
upstream in the cooling layer. For example, the Ha emissivity for fully developed radiative 
shocks increases approximately as Vj A . 

4.4. Observed Shock Velocities and Pre-Shock Densities 

Since the gas is relatively cool (T R3 10 4 K) in the photoionization zone, the optical 
lines have thermal widths of order lOA^ 1 / 2 km s~ x , where A is the element's atomic mass. 
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Macroscopic motion of the cooled layer will cause the observed line profiles of the optical lines 
to be significantly broadened. As the blast wave wraps around a protrusion, we expect to 
observe velocity components traveling both toward and away from us. The measured widths 
of the optical lines in the early HST /STIS spectrum suggests that the fastest radiative shock 
has a projected velocity ~ 250 km s" 1 (Michael et al. 2000). Moreover, the observed line 
ratios of Spot 1 indicate that lower velocity (< 135 km s" 1 ) shocks must also be present 
(cf. §5.2.1). Following the discussion in §4.1.1, we are not surprised that a range of shock 
velocities is required to explain the observations. 

Figure 14 shows the boundaries separating radiative and non-radiative shocks with t coo \ 
of 2 and 4 years overlayed on the hydrodynamically allowed regions of phase space for different 
values for the blast-wave pressure. In order for shocks with velocities as low as 135 km s _1 
to be radiative, the density of the obstacle p ;> 3 x 10 4 amu cm -3 has to be higher than 
most of the gas in the equatorial ring. A lower p is needed if the blast-wave pressure is 
lower than our nominal estimate. While faster shocks propagating into less dense gas may 
also be present in Spot 1, they do not contribute significantly to the observed optical and 
UV spectra because they are non-radiative. 

5. Analysis 

5.1. Nebular Analysis 

The HST /STIS spectrum of Spot 1 consists of many forbidden emission lines of various 
species in the UV and optical wavelengths (cf. Table 2 and 4). A standard nebular analysis 
on these data provides a good starting point because it indicates the typical values for the 
basic physical quantities that can be expected in the line-emitting region. We decide to focus 
our attention on the data taken around 1999 September (~ day 4570) when a complete set 
of spectrum is obtained from 1150 — 10270 A. For each ion we consider a five- level model 
and include all relevant atomic processes, such as collisional excitation and de-excitation, 
and spontaneous radiative transitions. The atomic data of Osterbrock (1989) are used for 
the analysis, except for N IV, where the results from Ramsbottom et al. (1994) are used. We 
construct a grid of line ratios as functions of electron number density n e and temperature 
T e . Figure 15 shows the contours for which ratios of lines of [O I], [O III], [N II], [S II], and 
[S III] agree with the observed values. 

The observed ratios of various forbidden lines require a relatively high electron density, 
n e ~ 10 6 cm" 3 , and temperature, T e ps 10 4 K. The upper limits for the line ratios of [N I] 
(5199 + 5201/3467 < 0.11) and [O II] (3726 + 3729/7320 + 7330 < 0.3) are consistent with 
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the derived high n e , as well as the N IV] (1483/1486) line ratio. Recalling that the typical 
number density of the ring is n ~ 10 4 cm -3 (Lundqvist & Fransson 1996), our result indicates 
that the gas in Spot 1 must have been compressed by a factor much greater than 4. This is 
consistent with our general picture (§4) that the UV and optical line emissions from Spot 1 
originate from a highly compressed gas behind a radiative shock. Our results also suggest the 
temperature and density stratification of the region where these emission lines are formed. 
This is also consistent with our general picture for radiative shocks where emission lines from 
different ionic species are created in different regions of the shock. 

As discussed earlier (§4.3), the optical emission lines from a radiative shock comes 
from a photoionization zone in the cooled post-shock gas. This photoionization picture is 
supported by the observed Balmer decrement in Spot 1, Ha : H/3 : H7 : = 3.63(±0.21) : 
1 : 0.48(±0.06) : 0.24(±0.04), which is generally consistent with both Case A and Case B 
values typically seen in photoionized gaseous nebulae (Osterbrock 1989). A complete analysis 
of the H and He lines will be presented in Pun et al. (2002, in preparation). 

5.2. UV Emission Line Modeling 

Emission lines from the higher ionization stages (e.g., N IV, C IV, O IV, N V) are formed 
in the cooling region of a shock, where the gas temperature and density are stratified and 
the ionization levels can be out of equilibrium. However, within each ionization stage, the 
population of all energy levels are in equilibrium. We model these lines with a radiative 
shock code that properly account for these effects. We compare our model results with the 
observations of six UV emission lines detected with the G140L grating (N V AA1238, 1243, 
Si IV A1393, 1403, O IV] A1400 multiplet, N IV] A1486, C IV AA1548, 1543, and He II A1640). 
We choose these lines because they have negligible contamination from the ring emission, 
and because they have been observed at three different epochs from day 3869 to day 4596. 
We decide not to model the Mg II AA2796, 2803 doublet, the strongest emission line in the 
Spot 1 spectrum, because of the large uncertainties in the observed flux (over a factor of 
2, cf. Table 4), caused by huge interstellar Mg II absorption along the LMC line of sight 
(cf. §3.3). Moreover, with the Mg II line flux measured only one epoch, we are not able 
to monitor its time evolution as for the other FUV lines. Similarly, we have not attempted 
models to fit both the UV and the optical emission lines because of the limited time and 
wavelength coverage of the optical data (cf. Table 2). 

For a given set of input shock-model parameters (shock velocity, pre-shock density, and 
abundances), we follow the evolution of the downstream gas until the shocked gas has cooled 
to 5000 K. We store the integrated flux for various emission lines at intermediate grid points 
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behind the shock as the basic vectors by which we compare with the observed data. We then 
vary the parameters in the models to obtain the best fit. To seek the best fit, we allow the 
abundances to vary from our standard "Ring" abundance described in §4.2. 

5.2.1. Single-Shock Models 

Using a grid of single-shock models with V s = 100 — 400 km s -1 , we were unable to 
reproduce the emission-line fluxes and widths observed from Spot 1 at any epoch. One 
of the main challenges for the single-shock model is to explain the relatively low ratios of 
N V A 1240 flux to the fluxes of lines from lower ionization stages (e.g., N IV] A 1486 and 
C IV A1550) in light of the widths (~ 250 km s" 1 ) of the observed lines. The abundance 
of N V increases rapidly for shock velocities exceeding th 135 km s _1 , while N IV and C IV 
are present at lower shock velocities (cf. Figure 13). As a result, the model line ratios 
N V A1240/C IV A 1550 « 10 and N V A1240/N IV] A1486 « 8 for shocks faster than 
135 km s^ 1 exceed the observed ratios of 0.6 — 1.7 and 2.5 — 5.6, respectively. The fact that 
this discrepancy exists for both line ratios indicates that this result does not depend on our 
choice of abundances. We therefore conclude that a single-shock model cannot reproduce 
the observed line emissions of Spot 1, as in our earlier studies (Michael et al. 2000). 

5.2.2. Two-Shock Models 

We are not too surprised that single-shock models fail to fit the observed emissions from 
Spot 1 because we expect there exists a range of shock velocities V s in the interaction (cf. 
§4.1.1). The exact distribution of V s depends on the geometry and density distribution in the 
obstacle, which are both unknown. Rather than attempting to explore the vast parameter 
space of possible shock velocity distributions, we fit the observed line emissions with models 
consisting of only two distinct shocks, each having a different velocity and surface area. 
We explore parameter space with shock velocities ranging from 100 to 400 km s _1 and find 
reasonable fits to the relative line strengths with a combination of one "slow" shock with 
V s = 135 km s _1 and one "fast" shock with V s = 150 - 400 km s -1 . As discussed below 
in §5.3, the observed line profiles suggest that the fastest radiative shocks have velocities 
V s «3 250 km s~\ We therefore assume this value for the fast shock component. 

The hydrodynamic simulations shown in Figure 8 of §4.1 illustrate two possible scenar- 
ios to account for the observed increase of the line fluxes with time. The crucial parameter 
distinguishing these scenarios is the ratio t s hock/t C ooh where t s h oc k is the "age" of the interac- 
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tion (i.e., the time since first impact), and t coo i is the characteristic radiative cooling time for 
the transmitted shocks. For the case where both the blast-wave pressure and the obstacle 
density are low, then t s h oc k ~ t coo i and the shocks in the obstacle would be in the process of 
developing their radiative layers. In this scenario the observed increase of the line fluxes is 
due primarily to development of new radiative layers as the shocks age. Alternatively, for 
the case where the density of the obstacle is high, then t s h oc k » tcooi and the shocks quickly 
become radiative. The observed increase of the line fluxes is then due to the increase in 
the surface area of the shocks as the blast wave overtakes more of the obstacle. We suspect 
that we are observing a combination of these two behaviors. Given our limited data set 
(6 lines in 3 epochs), we decide to fit only two limiting models which probably bracket the 
actual situation: evolution due solely to shock aging (Model 1), and evolution due solely to 
increasing shock areas (Model 2). 

For Model 1 we assume that the area of each shock remains constant over time and the 
only fitting parameters are the surface area and age of each of the two shocks. We assume 
that the faster shock is older because it represents the head-on shock that appears when the 
blast wave first encounters the obstacle, and that the slower shock is younger because it is 
driven into the sides of the obstacle by an oblique blast wave at a later time. For this model 
we need to assume a low pre-shock density n = 10 4 cm -3 so that t cool t shock 3 yrs 
for the fast shock. We also assume that the pre-shock density is constant throughout the 
obstacle for both models. The best-fit model for Model 1 has a Xr ~ 1-6 and its parameters 
are listed in Table 7. The model (squares) and observed (error bars) fluxes for all the 6 
emission lines are plotted in Figure 16. For the best-fit model, the surface area of the slower 
shock needs to be ~ 4 times greater than that of the faster shock to fit the low observed 
N V A1240/C IVA1550 flux ratio. In this model, the fast shock was 2.2 years old at the time 
of the first UV observation (day 3869), while the slower shock was only one week old. The 
age of the fast shock is only slightly smaller than the results of Lawrence et al. (2000a), who 
suggest that Spot 1 was at least 2.6 years old by the time of our first UV observation. This 
pleasant result is offset by the suspiciously young age of the slower shock. 

For Model 2 we assume that the pre-shock density is high enough so that all the shocks 
are fully radiative at all epochs. The fitting parameters more this model are the areas of the 
two shock surfaces at each of the three observation epochs. We find that the quality of the 
fits is rather insensitive to the assumed pre-shock density (no = 3.3 x 10 4 cm -3 is assumed). 
This is true because the 6 UV lines we are fitting are not subject to collisional supression 
(including N IV] A 1486, for which n crit > 3 x 10 9 cm -3 ), so their emission efficiencies, 77, do 
not depend on the pre-shock density (cf. §4.3). The luminosities of such lines are actually 
proportional to the product, n^A, of the pre-shock density and shock area. We present the 
best-fit paramters of the model in Table 7 and the model line fluxes in Figure 16 (circles). 
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Since Model 2 provides a better statistical fit (Xr ~ 1) to the observed line strengths than 
Model 1, we conclude that the increase of shock areas is the dominant cause of the increase 
of the fluxes. 

Figure 17 shows the time dependence of the fitted shock areas (actually the product 
n A) for Model 2. Similar to Model 1, we find that ~ 4 times more surface area must 
be covered by the slow shocks than by the fast shocks. While our best fit model results 
are statistically consistent with a constant ratio of slow to fast shock area (Ai35 : A 2 50 in 
Table 7), they suggest that this ratio may be decreasing with time. This behavior is not 
what we would expect from simple protrusion models shown in Figure 8. In such models 
the fast shocks are created first and only the area of slow shock interaction increases as the 
blast wave overtakes more of the obstacle. One possible explanation is that there may be 
many protrusions or clouds distributed in Spot 1 and therefore the net area of fast shocks 
can increase with time as more of them are encountered by the blast wave. Another possible 
explanation arises from the fact that Model 2 by construction does not allow for any aging 
of shocks. However, the real light curve of the Spot 1 emission is probably a combination of 
increased shock areas and additional emission from older, but newly cooled, radiative shocks. 
The addition of fast shock area at later times seen in Model 2 can therefore be attributed 
to shocks that were created before our first observation but added their emission later after 
they had cooled and become radiative. 

Our finding that more area is covered by the slow shocks than fast shocks suggests 
that the obstacle may have an elongated shape with an aspect ratio A s i ow /Af ast , where 
A s i ow and Af ast , are the shock areas for the slow and fast shocks, respectively. Assuming 
a simple, hemispherical geometry for the obstacle, we estimate that it has a typical scale 
length of d spot « [2(A s i ow + Af^A/ir] 1 / 2 . We find a best-fit size of the obstacle d spot « 
2 x 10 16 (3.3 x 10 4 cm^ 3 /n ) 1/2 cm at day 4606, the time of our last STIS observation. This 
is consistent with our finding (cf. §2) that Spot 1 remains spatially unresolved in our data. 

While the two-shock models successfully fit the time evolution of the UV emission lines, 
these models fail to account for the observed fluxes of the optical lines (e.g., Ha, H/3, [N II] 
AA6548, 6583, [O I] AA6300, 6364, and [O III] AA4959, 5007) emitted in the photoionization 
zone (§5.1). The observed optical fluxes are typically factors of 2 — 3 greater than those 
predicted by our models. We exclude emission from the precursor as a possible contributor 
to the extra emission because its emission would be narrow and center at zero velocity. 
However, we exclude that portion of the line profile when we fit the Spot 1 medium resolution 
G750M data (cf. Fig. 4). We discuss other possible explanations for this discrepancy in §6. 

Finally, we emphasize that the satisfactory fits of the UV line fluxes were possible only if 
we allow abundances to vary from our standard "Ring" abundance. The inferred abundances 
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are roughly consistent with the ring abundances inferred by Lundqvist & Fransson (1996), 
but differences as large as a factor of two are derived. 

5.3. Line Widths and Profiles 

In this section we discuss what can be learned from the widths and profiles of the 
observed emission lines. The data are sparse. We measure line widths of Spot 1 accurately 
for only a few optical lines (cf. §3.2). We also only have rough estimates of the widths of a 
few far UV lines (cf. §3.3.3). 

The line widths and profiles are dominated by the shock dynamics, and not thermal 
broadening. The line-emitting gas behind a radiative shock travels at approximately the 
shock speed and will therefore emit line radiation with a Doppler shift determined by the 
line-of-sight velocity of the shock. For emission lines originated from the photoionization zone 
(T « 10 4 K), e.g., the optical forbidden lines, the thermal broadening is ~ 10A- 1 / 2 km s~\ 
which is much smaller than the observed widths of the lines (~ 250 km s _1 ). The thermal 
broadening for emission lines formed in the cooling region, such as N V A1240, while larger 
than that for the optical forbidden lines, is still much smaller than the dominating Doppler 
broadening effects of the macroscopic motions of the shock. Although it is poorly measured, 
the width of N V A1240 line (Table 4) appears to be higher than that of all the other optical 
and UV emission lines, with the possible exception of Ha. This result can also be explained 
by Figure 13, which shows that the N V A1240 emission can only be produced in shocks 
faster than 135 km s _1 , while all the other lines can be produced by even slower shocks. 

The fact that the measured widths of the [O I], [N II], and [S II] lines are smaller 
than Ha (Table 3) indicates that the faster shocks produce these lines with lower efficiencies 
(compared to Ha) than the slower shocks. This is caused by collisional supression of the 
forbidden lines, which increases with increasing compression ratio of the photoionized zone, 
as shown previously in Figure 13. 

The line profile shapes depend on the geometry of the shock surface and its orientation 
relative to the observer. The shocks having the highest projected velocities toward and away 
from the observer produce the blue and red wings of the lines, respectively. The Doppler 
shifts of the wings of the line profile therefore provide a lower limit to the velocity of the 
fastest radiative shocks present while both fast transverse shocks and slow shocks contribute 
to the cores of the line profiles. 

We construct simulated Ha line profiles based on simple geometries for the protrusion 
(e.g., Figure 8) in order to derive constraints on the hydrodynamics model of the interaction. 
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The 1999 August 30 (day 4568) G750M 0"1 line profile of Ha is chosen because it is well- 
observed and is not affected by collisional supression. Figure 18 shows our model shock 
surface geometries. For each model we assume a distribution of shock velocities (from 100 to 
250 km s _1 ) normal to the surface of the protrusion. We normalize each velocity distribution 
by requiring that the relative areas of fast (V s > 135 km s _1 ) and slow (V s < 135 km s _1 ) 
shocks in each model is ~ 4, as determined from the best-fit two-shock models (§5.2.2). 
We also assume that the axis of symmetry of the protrusion points radially inward from 
the circumstellar ring at the location of Spot 1 for all models. We calculate line profiles by 
convolving the surface emissivity (oc V^ A , cf. §4.3) of the radiative shocks with the fractional 
areas covered by shocks having the appropriate line-of-sight velocities. We then convolve the 
model profile with the STIS LSF (Leitherer et al. 2000). At the Ha wavelength, the STIS 
LSF for the G750M O'.'l aperture can be well approximated with a Gaussian of FWHM 1.5 
pixels. The simulated profiles are finally normalized for comparison with the observed Ha 
profile. 

Figure 19 shows the resulting Ha line profiles (solid) from the geometry models of 
Figure 18. For each case the day 4568 Ha profile (squares) is also plotted for comparison. 
The observed line width indicates that radiative shocks are present with projected velocities 
« 250 km s _1 , which is a good indication of the velocities of the fastest radiative shocks 
present. The line profile formed from an ellipsoidal model surface (model A) has multiple 
peaks and steep edges that are not apparent in the observed profiles. The peaks are due to 
caustics in velocity space that occur where the shock surface lies tangent to the plane of the 
sky. These features are less pronounced in the line profiles produced by a spherical model 
surface (model C), where tangential surfaces comprise of a lower fraction of the total surface 
areas. 

For each case we also show a model in which no radiation comes from the tip of the 
protrusion (models B and D, for elliptical and spherical surfaces, respectively). These models 
are motivated by the observations from the Chandra X-Ray Observatory (Burrows et al. 2000; 
Park et al. 2001), which show a bright spot in the X-ray image of SN 1987A coincident with 
Spot 1. One possible interpretation is that the X-rays come from a faster (V s > 600 km s _1 ) 
non-radiative shock at the tip of the protrusion responsible for Spot 1. We notice that this 
picture is also consistent with one where there is a density gradient across the ring. 

While the model line profiles have similar widths as the observed Ha line, they do not 
provide a good fit in the high velocity portion (;> 100 km s _1 and <^ —150 km s _1 ) of the 
line profile. As we shall discuss in §6, radiative shocks are subject to instabilities that will 
introduce "turbulence" in the radiative layer. We simulate this effect of this turbulence qual- 
itatively by convolving the derived model line profiles with Gaussians of FWHM 100 km s _1 . 
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The resulting line profiles (dashed lines in Figure 19) provide substantially better matches 
to the observed Ha line. 



6. Discussion 

While the idealized radiative shock models described in §4 and §5 provide a plausible 
framework for interpreting the spectrum of Spot 1, they do not account for all the details of 
the observed line ratios and profiles. In this section, we consider what the departures of the 
actual spectra from the model spectra may tell us about how the actual physical situation 
differs from our models. 

First, we discuss the "turbulent" broadening that we introduce in §5.3 to make the model 
line profile resemble the observed one at high velocity. This turbulence is almost certainly due 
in part to the cooling instabilities that exist behind radiative shocks with V s j> 150 km s _1 
(Chevalier & Imamura 1982; Imamura et al. 1984; Innes et al. 1987; Strickland & Blondin 
1995; Walder & Folini 1998). To illustrate the effects of these instabilities, we construct a one- 
dimensional VH-1 simulation of a blast wave driving a radiative shock into gas of constant 
upstream density, no = 3.3 x 10 4 cm -3 and the results are shown in Figure 20. In the figure, 
the upstream gas flows to the left with a velocity V in fi ow = 220 km s _1 while the radiative 
gas, shown as a dark band, propagates to the right with a mean velocity ~ 250 km s -1 . The 
dark gray parabolas represent the positions of the actual shock front, which moves ahead 
of the radiative layer and then collapses back to the radiative layer over time. The velocity 
of the shock front relative to the upstream gas varies from ~ 100 to 300 km s _1 . This 
process repeats itself on a time scale comparable to the cooling time of the shocked gas. In 
our simulation, the forward shock never stabilizes while the oscillation of the velocity of the 
radiative layer decreases as it gains mass with time. 

The instabilities that exist behind radiative shock become more dramatic in the two- 
dimensional simulations, as described in Walder & Folini (1998). We construct a 2-D VH-1 
shock simulation in Figure 21 to illustrate such effects. We seed the instability with a 5% 
sinusoidal density perturbation of the unshocked gas in a direction perpendicular to the shock 
front, with the purturbation wavelength comparable to the cooling length of the transmitted 
shock. As the shock enters the dense gas, a Richtmyer-Meshkov instability develops rapidly 
at the contact discontinuity. After several cooling times (and forward shock oscillations), 
the radiative layer is shredded into a chaotic flow of tiny blobs by a combination of thermal 
and Rayleigh- Taylor instabilities. 

If the blast wave strikes the obstacle at oblique incidence, the flow of shocked gas 
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parallel to the contact surface will cause strong Kelvin-Helmholtz instabilities to shred the 
surface even further and will introduce a component of velocity parallel to this surface in 
the radiative layer. We can see these instabilities beginning to develop in the last frames of 
Figure 8. However, we have not attempted to model their consequences in our simulations 
of the line profiles. 

Since the HST only allows us to resolve structures with dimensions J> 10 17 cm at 
the distance of SN 1987A, we have no direct knowledge of the actual shape or the density 
distribution of the protrusion gas responsible for the Spot 1 emission. Even in the case where 
it has a simple shape such as that shown in Figure 8, the instabilities discussed above will 
introduce turbulent structures at scale lengths ~ 10 13 cm, which we cannot hope to resolve 
with current technology. 

Moreover, it is unclear whether the assumption that the protrusion gas has a simple 
shape with a smooth surface is justified. While we know that the circumstellar ring has 
a range of densities (Lundqvist & Fransson 1996), we do not know the scale lengths upon 
which these inhomogeneities may be distributed. If the blast wave encounters a "lumpy" 
medium, a complex flow will develop where the radiating blobs are hammered both by the 
reflected shocks and by the blast wave itself. Clearly, the current observations allow for many 
plausible scenarios for the hydrodynamics of the radiating gas for which quantitative models 
are beyond the scope of this paper. We are not surprised that the simple hydrodynamic 
models we have considered in §5.3 fail to reproduce the observed line profiles. 

Magnetic fields may also affect the shock dynamics. For the range of densities and shock 
velocities expected in Spot 1, we estimate that a magnetic field with strength Bq j> 400 /xG 
entrained in the upstream gas would be sufficient to suppress compression in radiative shocks 
and stabilize the cooling instabilities (Innes 1992). This will result in a photoionization zone 
of lower density. 

In our idealized shock models of §4 and §5, we have only considered radiative shocks 
in which the pressure driving the shocks is constant. However, as Borkowski et al. (1997b) 
and Luo et al. (1994) have shown, the pressure behind an actual blast wave encountering 
an obstacle have a complex time dependency. As the reflected blast wave travels away from 
the surface of the obstacle the flow diverges and the driving pressure decreases. Later, the 
pressure driving the radiative shock into the obstacle may suddenly rise again owing to 
"echo" shocks produced when the reflected shock strikes the contact discontinuity and the 
reverse shock. 

The most troubling failure of our models is their inability to account for the observed 
ratio of fluxes of optical to UV lines. As we have described, the radiative shock models 
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consistently predict optical fluxes that are factors R3 3 lower than expected from models that 
are normalized to give the observed UV line fluxes. 

A simple way to account for the discrepancy between the observed and model opti- 
cal/UV line ratios is to assume that the extinction of Spot 1 is greater than that previously 
determined for SN 1987A. If we applied a reddening of E(B — V) « 0.23, the corrected 
UV flux would increase by a factor 3, sufficient to remove the discrepancy. Consider- 
ing that Spot 1 is located on the far side of the ring from the observer, a column density 
n£ ~ 10 21 cm" 2 of gas with a 30 Doradus dust/gas ratio in the foreground of Spot 1 (perhaps 
the ring material itself) could account for such additional reddening. It is also possible that 
the shock itself may destroy graphite dust grains and elevate the ratio of C/N abundances in 
Spot 1 (compared to the unshocked ring) to account in part for the relatively low observed 
ratio of N V A1240/C IV A1550. Shocks faster than 200 km s" 1 can destroy graphite 
grains in a timescale ^ 2 yrs (Tielens et al. 1994). 

The revised value of E(B — V) is 0.04 (and 2a) greater than the measurement of 0.19 ± 
0.02 by Scuderi et al. (1996) on the nearby Star 2, and 0.07 greater than the value we adopted 
for Table 4. However, modeling of the UV and optical narrow emission lines for the inner 
circumstellar ring by Lundqvist & Fransson (1996) does not support such large E(B — V) 
value. Therefore we do not regard it to be highly probable. 

A second way to account for the optical/UV line ratio is the presence of an extra source of 
ionization in the photoionization zone that is roughly three times stronger than the radiation 
from the shocked gas itself. In the one-dimensional radiative shock models described in §4, 
the only source of ionizing radiation for the photoionized zone is the shocked gas immediately 
upstream from that layer. One possible candidate for the extra ionization is the faster shocks 
which have not had time to develop their own radiative layers. However, we are not confident 
that such models will work. The problem is that the presumed extra ionizing radiation 
must also illuminate the cooling zone where the UV emission lines are formed, and this 
illumination will elevate the abundance of N V in the cooling layer. As a result, shocks with 
velocity V s <^ 135 km s -1 will produce extra N V A1240 emission, causing the predicted N V 
A 1240 flux to again exceed the observed value. A more promising candidate for the extra 
source of ionization is the nonthermal radiation from particles accelerated in the shock itself. 
Detailed modeling is required to demonstrate its feasibility. 

A third possible explanation for the failure of the models to account for the ratios 
of optical to UV emission lines is the inadequacy of the stationary radiative shock models 
themselves. Our findings are similar to that from the modeling of the spectra of Herbig-Haro 
objects with simple plane-shock and bow-shock models, where high-ionization lines, such as 
C IV and N V, are underproduced (Bohm & Goodson 1997). As we have shown in Figure 21, 



-32 - 



radiative shocks are inherently unstable and the ratio of optical to UV emission lines from 
such unstable shocks may be different from that in a stationary flow, even in a time-averaged 
or ensemble-averaged sense. For example, the temperature gradients in the unstable flow 
are more pervasive than in a stationary flow, in which case thermal conduction may provide 
an additional energy source. 

Future observations with the HST may be very helpful in resolving some of the issues 
that we have raised here. Valuable clues to resolving the optical/UV line ratio puzzle may 
come from observations of the same ratio in the other hot spots. For example, if the ratio is 
affected by extinction by dust in the immediate foreground of Spot f , it might have different 
values in the other spots, some of which lie on the near side of the circumstellar ring. 

We can also benefit from observations of the UV emission lines at a higher spectral 
resolution. First, we can better estimate the amount of flux reduction caused by the LMC 
interstellar absorption (cf. §3.3.3) by obtaining the line profile shapes of emissions not 
affected by this absorption, such as N V A1240, and N IV AAI483, 1486. This will significantly 
reduce the error estimates of the measured fluxes affected by this absorption, especially for 
Mg II A2800, where the flux correction due to LMC absorption applied is both the highest 
and the most uncertain. Second, to reconcile the observed low flux ratio of N V A1240/C IV 
A1550, we have suggested a two-shock model in which the higher- velocity radiative shocks 
make up a relatively small fraction of the shock area compared to the lower-velocity shocks. 
This model implies that the ratio of N V A1240/C IV A1550 line emissivity should increase 
in the wings of the lines, where only the fastest shocks contribute. On the other hand, we 
have also suggested that destruction of graphite grains in the faster shocks may elevate the 
abundance of carbon in the gas phase and lead to a decrease in the ratio of N V A1240/C IV 
A1550 in the wings of the lines. However, as shown in Figure 6, the observed profile of C IV 
A1550 is severely modified by interstellar absorption, introducing large uncertainties in the 
N/C abundance ratio. We may get a clear look at the red wing of C IV A1550 by observing 
the spots on the far side of the ring (e.g., Spots 3, 4, and 5), where the shock emission would 
be mostly redshifted to wavelengths where the interstellar absorption is not so severe. When 
Spot 1 becomes bright enough to permit observations of the line profile of N IV] A 1486, the 
flux ratio of N IV] A1486/N V A1240 as a function of position within the line profile will give 
us a diagnostic of shock excitation versus velocity that is independent of relative abundances 
and not subject to interstellar absorption. 



7. Summary 

In conclusion we summarize the main results of this paper: 
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• We describe the HST/STIS UV and optical spectrum of Spot 1 up to day 4606 (1999 
October 7). We measure and tabulate fluxes, line widths, and line centroids for all 
lines detected in the spectrum. 

• We determine correction factors to account for the substantial interstellar line absorp- 
tion of several lines (C II A1335 multiplet, Si IV AA1394, 1403, C IV AA1548, 1551, 
and Mg II AA2796, 2803). However, these factors are poorly constrained because of 
uncertainties in the line widths and centroids. 

• The observed emission is caused by radiative shocks which develop when the supernova 
blast wave strikes dense gas protruding inward from the equatorial ring. A nebular 
analysis of several line ratios confirms that many of the emission lines are formed in a 
region of very high density gas (n e ~ 10 6 cm -3 ), which has been compressed by factors 
> 100 due to radiative cooling downstream from the shock. 

• The observed line widths indicate that radiative shocks with velocities as high as 
250 km s -1 are present, while the line ratios indicate that much of the radiation must 
come from slower shocks (V s < 135). The inferred range of shock velocities is a natural 
consequence of a model in which a blast wave overtakes a dense obstacle. 

• Hydrodynamic arguments show that such slow shock velocities will be present only if 
the density in the unshocked obstacle is sufficiently high (n ~ 3 x 10 4 cm" 3 ). This 
density lies at the high end of the range of densities observed for the ring by Lundqvist 
& Fransson (1996). The fact that shocks with velocity V s = 250 km s -1 must have 
become radiative after 3 years sets a lower density limit, n ^ 10 4 cm -3 . 

• Faster shocks (V s f=a 250 - 1000 km s -1 ) may also be present in Spot 1, as would be 
the case for shock interaction with the lower-density gas (n < 10 4 cm -3 ) present in 
the ring (Lundqvist & Fransson 1996). These shocks would be invisible in the optical 
and UV emission spectra though because they would not yet have developed radiative 
layers. However, they may contribute to the enhancements of X-ray emission seen 
near the hot spots by the Chandra X-ray Observatory (Burrows et al. 2000; Park et 
al. 2001). 

• We interpret a subset of the observed emission lines (6 UV lines observed at 3 epochs) 
with a model consisting of two one- dimensional stationary radiative shocks with dif- 
ferent velocities (V s = 135 and 250 km s -1 ). In our favored model, the lines brighten 
because the area of the obstacle overtaken by the blast increases with time. The area 
covered by the slower shocks must be « 4 times that covered by the faster shocks. To 
account for the observed line strengths, the obstacle must have a characteristic scale 
length of 2 x 10 16 cm. 
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• Our models underestimate the ratio of observed optical/UV line fluxes by a factor of 
~ 3. This may indicate the presence of an additional source of ionizing flux in the 
photoionization zone, or that the assumption of steady state shocks is too limiting. 

• The observed line profiles are in qualitative agreement with simulated profiles created 
from simple geometries for the obstacle. The nearly Gaussian observed profiles indicate 
the presence of chaotic flows in the photoionized region, caused by violent thermal and 
shear instabilities in the shocked gas. 

• Continued UV and optical observations of Spot 1 are necessary to further our under- 
standing of the physics of the radiative shocks present there. Resolving the profile of 
the N V A1240 and other UV lines will provide critical tests of the radiative shock 
model. 

• Observations and modeling of the emission-line spectra of the other spots now appear- 
ing in the ring will provide vital clues to resolving the uncertainties of the present 
model for Spot 1. Observations with HST/STIS are required to spatially resolve their 
individual spectra. 
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Table I. HST/STIS observations of Spot 1. 



Grating 






Wavelength 


Resolution 




Days after 


Exposure 


setting 


Slit 


Range (A) 


(A) 


Date 


outburst 


time (sec) 


G140L (1425) 


52" 


x 0'.'5 


1140- 1700 


1.0 


1997 September 27 


3869.3 


11222 


G140L (1425) 


52" 


x 0"5 


1140 - 1700 


1.0 


1999 February 27 


4387.6 


14350 


G140L (1425) 


52" 


x 0"2 


1140- 1700 


1.0 


1999 October 7 


4606.1 


10478 


G230L (2376) 


52" 


x 0"2 


1570 - 3180 


3.0 


1999 September 17 


4586.5 


10125 


G430L (4300) 


52" 


x 0"2 


2900 - 5700 


4.0 


1999 September 3 


4571.9 


7583 


G750L (7751) 


52" 


x 0"5 


5240 - 10270 


8.0 


1999 February 21 


4380.9 


10500 


G750L (7751) 


52" 


x 0"2 


5240 - 10270 


8.0 


1999 September 18 


4387.4 


7583 


G750M (6581) 


52" 


x 0"2 


6295 - 6867 


1.0 


1998 March 7 


4030.0 


8056 


G750M (6581) 


52" 


x 0"1 


6295 - 6867 


1.0 


1999 August 30 


4368.0 


3 x 7804 a 



a Not all observations centered on Spot 1 (refer to text). 
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Table 2. HST/ST1S dereddened optical emission-line fluxes from Spot 1 





Central 


1998 Mar 7 


1999 Feb 21 


1999 Sep l a 


1999 Sep 18 




Emission 


Wavelength 


(day 4030.0) 


(day 4380.9) 


(day 4570.0) 


(day 4587.4) 


Extinction 




(A) 


a b 

flux D 


n h 

flux D 


flux D 


flux" 


correction 


rivr-, Vl 








O Q -L 1 P. 




9 HQ 


[N Tl 
I J-J 


3466.5 






7.2 ± 1.9 




2.08 


ro III 


^79^ fl 4- ^798 8 






< 6.0 




2.02 


Hi 


3770.6 






1.2 ± 1.0 




2.01 


H0 


3797.9 






1.6 ± 1.2 




2.00 


Hr? 


3835.4 






4.4 ±1.3 




2.00 


[Nc III] 


3868.7 






13.3 ± 2.7 




1.99 


He I + Hf 


3888 7 4- 3889 1 






18 + 2 8 




1.98 


Ca II + [Fc III C8F1 


3933 7 4- 3931 4 4- 3932 7 






4.6 ± 1.0 




1.97 


He I + Ca II + [Ne III] 


3964.7 + 3968.5 + 3968.7 






12.5 ± 2.6 d 




1.96 


He 


3970.1 










1.96 


He I 


4026.2 






2.3 ± 1.0 




1.95 


fS III 


4068.6 






21.8 ± 2.4 




1.93 


IS III 
L° - LJ -J 


4076.4 






10.5 ± 1.3 




1.93 


H5 


4101.7 






12.9 ± 1.8 




1.92 


[Fe II] 


4244.0 






2.6 ± 0.8 




1.88 


X17 


4340.5 






24 q + 2 5 




1.85 


[F e III (21F) 


4358.4 






4.9 ± 2.5 d 




1.84 


fO IIIl 

J 


4363.2 






4.9 ± 2.5 d 




1.84 


He I 


4387.9 






1.2 ± 0.5 




1.84 


[Fe III (6F) 

L J V w / 


4416.3 






3.2 ± 0.7 




1.82 


[Fc II] (7F + 6F) 


4452.1 + 4458.0 






2.0 ± 0.7 




1.81 


He I + [Fc II] (6F) 


4471.5 + 4470.3 






4.6 ± 0.7 




1.81 


Mg I] 


4571.1 






4.1 ± 1.0 




1.78 


[Fe II] (4F + 5F) 


4664.5,4665.7 + 4665.0 






2.0 ± 0.8 




1.75 


He II + [Fc II] (5F) 


4685.7 + 4687.6 






4.0 ± 0.7 




1.75 


[Fe II] (20F) 


4814.6 






1.7 ± 0.5 




1.71 


H/3 


4861.3 






49.6 ± 2.7 




1.70 


[Fc II] (3F + 4F) 


4889.7 






1.7 ± 0.7 




1.70 


[Fe II] (20F) 


4905.4 






1.4 ± 0.5 




1.69 


He I 


4921.9 






0.6 ± 0.3 




1.69 


[O III] 

J 


4958.9 






7.7 ± 1.4 




1.68 


[0 III] 


5006.8 






24.4 ±1.1 




1.67 


He I 


5015.7 






4.5 ±0.7 




1.67 


[Fe II] (18F + 19F) 


5158.0 + 5158.8 






4.6 ±0.6 




1.64 


[N I] 


5197.9 + 5200.3 






< 0.8 




1.63 


[Fc II] (19F + 18F) 


5261.6 + 5268.9,5273.4 






4.8 ±0.7 




1.62 


[Fe II] (19F) 


5333.7 






1.0 ±0.3 




1.60 


[Fe II] (19F) 


5376.5 




0.6 ±0.4 


1.4 ±0.3 


0.6 ± 1.0 


1.60 


He II + [Fc II] (17F) 


5411.5 + 5412.6 




1.0 ±0.3 


1.5 ±0.4 


1.4 ±0.6 


1.59 


[Fe II] (34F) 


5527.3 






1.1 ±0.3 


1.4 ±0.5 


1.57 


[O I] 


5577.4 




1.0 ±0.5 


1.9 ±0.4 


2.3 ±0.6 


1.57 


[N II] 


5754.6 




16.9 ±2.0 




20.8 ± 1.4 


1.54 


He I 


5875.7 




9.3 ±2.1 




14.4 ± 1.1 


1.53 


[O I] 


6300.3 


13.2 ± 1.0 


18.7 ± 1.5 


27.1 ± 1.2 


25.9 ± 1.9 


1.48 


[s III] 


6312.1 








0.6 ±0.3 


1.48 


[O I] 


6363.8 


4.9 ±0.6 


7.4 ±0.5 


8.1 ±0.9 


8.0 ± 1.3 


1.47 
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Table 2 — Continued 





Central 


1998 Mar 7 


1999 Feb 21 


1999 Sep l a 


1999 Sep 18 




Emission 


Wavelength 


(day 4030.0) 


(day 4380.9) 


(day 4570.0) 


(day 4587.4) 


Extinction 




(A) 


fiux b 


flux b 


fiux b 


fiux b 


correction 


[N II] 


6548.1 


8.2 ± 1.1 




9.5 ± 1.0 




1.45 


Ha 


6562.8 


67.9 ±1.8 




179.8 ± 3.9 




1.45 


[N II] 


6583.5 


19.0 ± 1.3 




31.1 ± 1.7 




1.45 


He I 


6678.2 


1.2 ±0.3 


1.6 ±0.2 


2.2 ± 0.5 


3.5 ±0.6 


1.44 


[s II] 


6716.5 


1.3 ± 1.3 




1.3 ± 0.8 




1.44 


[s II] 


6730.8 


1.3 ±0.5 




2.2 ± 0.7 




1.44 


He I 


7065.2 




4.3 ±0.9 




6.3 ±0.7 


1.40 


[Ar III] 


7135.9 




2.6 ±0.4 




2.5 ±0.7 


1.40 


[Fe II] (14F) 


7155.2 




3.0 ±0.6 




4.3 ±0.8 


1.39 


[Fe II] (14F) 


7172.0 




3.7 ±0.7 




1.2 ±0.5 


1.39 


[Ca II] 


7291.5 




4.1 ±0.7 




6.2 ±0.9 


1.38 


[Ca II] + [O II] 


7323.9 + 7320,7330 




20.9 ±4.0 




25.2 ± 1.3 


1.38 


[Ni II] (2F) 


7377.9 




0.4 ±0.5 




2.1±0.9 C 


1.37 


[Fe II] (14F) 


7388.2 




0.3 ±0.4 




2.1±0.9 C 


1.37 


[Ni II] (2F) 


7411.6 




0.4 ±0.4 




WEAK 


1.37 


[Fe II] (14F) 


7452.6 




0.4 ±0.2 




1.1 ±0.5 


1.36 


[Fe II] (13F) 


8617.0 




1.1 ±0.4 




4.4 ±0.9 


1.27 


[S III] 


9068.6 




0.8 ±0.5 




1.4 ±0.9 


1.24 


[S III] 


9530.6 




4.8 ±1.7 




11.2 ± 1.6 


1.22 



Combining the G750M data taken on August 30, 1999 (day 4568.0) and the G430L data taken on September 3, 1999 
(day 4571.9). 

b Flux in units of 10 -16 ergs cm~ 2 s _1 . 

C E(B — V) = 0.16, and the extinction correction law of Cardelli et al. (1989), with Ry = 3.1. 

d Emission blended in the G430L data, separate component determined from Jan 2000 data (Pun et al. 2002, in prepa- 
ration). 

°Emissions blended in the G750L data. 
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Table 3. Peak emission velocity and FWHM width of Spot 1 emissions 



Emission 


1998 Mar 7 (day 4030.0) 


1999 Aug 


30 (day 4568.0) 


Line 


Width a 


Peak velocity 3. 


Width a 


Peak velocity 3 


[0 I] A6300 


154 ± 7 


-13 ± 18 


172 ± 5 


-34 ± 5 


[0 I] A6363 


179 ± 14 


-27 ± 20 


182 ± 11 




[N II] A6548 


165 ± 11 


-37 ± 24 


230 ± 16 


-11 ± 4 


Ha A6563 


253 ± 4 


-29 ± 3 


225 ± 3 


-35 ± 2 


[N II] A6584 


204 ± 7 


-17 ± 10 


215 ± 7 


-19 ± 2 


He I A6678 


177 ± 38 




183 ± 27 




[S II] A6717 


113 ± 40 




84 ± 24 




[S II] A6731 


149 ± 48 




134 ± 20 




a Velocity in 


units of km 


s" 1 . 







-44- 



Table 4. HST/STIS dereddened UV emission-line fluxes and line widths of Spot 1 



Central Flux 1997 Sep 27 (day 3869.3) 1999 Feb 27 (day 4387.6) 1999 Sep 27 (day 4596.2) Extinction 



Emission 


Wavelength (A) a 


correction 15 


Flux c 


Width d 


Flux c 


Width d 


Flux c 


Width d 


correction 


nV 


1238.8 + 1242.8 


2 6+ 1 ' 5 


15.0 ± 5.3 


270 +220 


51.5 ±6.4 


380111° 


79.6 ± 5.6 


33011- 


5.93 


ClI 


1335 multiplet 










7 o + 5.4 
<■•} _ 4.2 




4.96 


Si IV 


1393.8 + 1402.8 


2 0+°' 8 


r: o + 3.2 
0.0 _ 2 g 




14.9 ± ti 


130 +150 


21-0 tVs 


160l 180 


4.57 


oIV] 


1400 multiplet 




6.1 ± 1.7 




11.0 ±2.4 




11.1 ± 2.9 




4.53 


NlV] 


1483.3 




5.9 ± 1.5 




6.5 ± 1.3 


110+ 130 


9.2 ± 1.4 


200+ 100 


4.16 


NlV] 


1486.5 




6.1 ± 1.5 




11.5 ± 1.9 


110+ i30 


14.2 ± 1.8 


200 +ioo 


4.15 


cIV 


1548.2 + 1550.8 




24 1 + 13 ' 3 
"^.1 _ g.o 




o 7 o +18.9 
-11.5 


150 +ioo 


45 7 +23.4 
40 - ' -14.4 


180 +i50 


3.95 


[Ne IV] 


1601.5 + 1601.7 












3.9 ±2.4 




3.81 


He II 


1640.4 




13.1 ±3.4 




36.7 ±4.8 


100 +80 


64.3 ± 8.9 


130 +80 


3.73 


OlII 


1660.8 + 1666.2 












13.7 ±2.8 




3.69 


ClII] 


1908.7 












19.7 ±5.2 




3.56 


Nil] 


2139.0 + 2142.8 












18.5 ± 3.2 




4.02 


ClI 


2325 multiplet 












33.4 ±2.2 




3.43 


Si II 


2334.6 












6.9 ± 1.5 




3.39 


[oil] 


2470.3 
2736.7 f 
2746.3* 












10.7 ± 1.5 
2.3 ±0.8 
3.0 ±0.8 




2.93 
2.51 
2.50 


Mg II 


2795.5 + 2802.7 












230 tf 3 l 




2.45 



a Wavelength in A vac for A < 2000 A, A air for A > 2000 A. 

b Corrected for LMC line absorption assuming Spot 1 FWHM velocity V FW HM = 1501^° km s" 1 and peak velocity V = -30 ± 15 km s" 1 . 
c Flux in units of 10 -16 ergs cm~ 2 s _1 . 
d Emission-line widths in units of km s -1 . 

C E(B — V)lmc = 0.06 and E(B — y)Galactic = 0.10 are assumed, along with Ry = 3.1. For the LMC component, the 30 Doradus extinctior 
function of Fitzpatrick (1986) is used; for the Galactic component, the Seaton (1979) correction law is used. 

'Uncertain identification, refer to text. 
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Table 5. Two-shock model results 



Model 


Model 1 


Model 2 


Parameters 3, 


no = 10 4 cm~ 3 


n = 3.3 x 10 4 cm- 3 


T"i35 (years) 


0.03 ± 0.01 




7250 (years) 


2.21 ± 0.01 




^4l35 : ^250 


3.6 ± 0.7 




^135 : ^250(1) 




5.9 ± 3.3 


^135 : ^250 (2) 




3.3 ± 1.2 


^135 : ^250 (3) 




2.8 ± 0.9 


He 


1 


1 


C 


1.85 ± 0.59 


2.06 ± 0.66 


N 


0.55 ± 0.08 


0.54 ± 0.08 





0.83 ± 0.16 


0.75 ± 0.14 


Si 


1.15 ± 0.37 


0.78 ± 0.24 



a The model parameters from top to bottom: (Model 1) 
age of each shock at the first observation epoch, ratio of 
shock surface (Model 2) ratio of shock surface areas 

for the three observations (day 3869, day 4388, day 4596); 
(both Models 1 and 2) abundances adjustment factor with 
respect to the "Ring" abundances (cf. §4.2). The He 
abundance is kept fixed for both models. 
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SN Center 
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Fig. 1. — Schematic representation of the double-shock structure of SNR 1987A. 
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Fig. 2.— (a) Slit position of 1999 February 27 G140L 075 observation plotted on the 1999 
April 21 WFPC2 Ha image of SN 1987A. The image displayed is a 4"56 x 4"56 square, (b) 
Same as (a), for the 1999 October 7 G140L 0'.'2 observation, (c) Same as (a), for the 1999 
September 17 G230L observation (solid) and for the 1999 September 3 G430L and 1999 
September 18 G750L 0"2 observations (dotted), (d) Same as (a), for the 1999 August 30 
G750M O'.'l observations. All three O'/l slit positions are plotted. 
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Fig. 3.— (a) Subsection of the 1999 October 7 STIS G140L (1425) spectral image showing 
the emission from Spot 1. The height of the image displayed is 6" along the slit, (b) Same 
as (a), for the 1999 September 17 STIS G230L (2376) data, (c) Same as (a), for the 1999 
September 3 STIS G430L (4300) data, (d) Same as (a), for the 1999 September 18 STIS 
G750L (7751) data, (e) Same as (a), for the 1999 August 30 STIS G750M (6581) data. 
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Fig. 4.— The profile of the [N II] A6583 Spot 1 emission obtained with the G750M (6581) 
grating in 1999 August. The low-velocity region (shaded) was dominated by emission from 
the circumstellar ring and was excluded from the fit. The best- fit Gaussian (dotted) to the 
Spot 1 emission is shown. The velocity scale has been corrected for the heliocentric velocity 
of the LMC (285 km s _1 ). The residual to the best-fit Gaussian is also shown. 
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Fig. 5. — (a) The interstellar absorption correction factors for emission lines C II AA1335,1336 
(dot-dash), Si IV AA1394, 1403 (dash), C IV AA1548, 1551 (dot), and Mg II AA2796, 2803 
(solid). The peak velocity, Vq, of the model Spot 1 emission is assumed to be —30 km s _1 . 
(b) The correction factor for C IV AA1548, 1551 doublet is plotted against the FWHM 
velocity of the input model Spot 1 Gaussian profile with peak velocities Vq — km s _1 
(solid), —25 km s _1 (dots), —50 km s" 1 (dash), —75 km s _1 (dot-dash), and —100 km s _1 
( dot- dot- dot- dash) . 
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Fig. 6. — (a) The IVE high dispersion spectrum of the C IV AA1548, 1551 absorption lines 
(solid), overlayed with a model Gaussian Spot 1 line profile of the C IV emission line (dots). 
The FWHM width of the model profile is 150 km s" 1 and it peaks at —30 km s" 1 with line 
ratio 7(1548) //(1551) = 2:1. (b) The model C IV Spot 1 spectrum after the interstellar line 
absorption, (c) The LSF of STIS/MAMA G140L grating for the 072 aperture at 1500 A. (d) 
The model "observed" spectrum after convolving the model profile after absorption with the 
line spread function (thick). Example of a Monte-Carlo generated line profile of the C IV 
doublet, normalized to the integrated counts of the 1999 October data, is also shown (solid), 
along with the best-fit double Gaussian function to this profile (dash). In this example, the 
best-fit Vfwhm = 281 km s -1 . 
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(d) C IV 154X 1551 
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Fig. 7. — Monte Carlo results of the relation between the intrinsic Spot 1 emission width and 
the measured widths from Gaussian fittings for the 1999 October G140L 0'/2 observations in 
(a) N V AA1239, 1243; (b) Si IV AA1394, 1403; (c) N IV] AA1483, 1486; (d) C IV AA1548, 
1551; and (e) He II A1640. 
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Scenario 1 Scenario 2 




Fig. 8. — 2D cylindrically symmetric hydrodynamical simulations (800 x 600 zones) of a blast 
wave overtaking a dense obstacle with darker shade representing higher density. The left pan- 
els (Scenario 1) show the development of a shock (V b = 2, 000 km s^ 1 , p Hn = 150 amu cm -3 , 
d sp ot = 4 x 10 16 cm, p = 10 4 amu cm -3 ) where the cooling time is comparable to the time 
for the fast shock to cross the obstacle, i.e., t coo i t cross . Shocks that have undergone 
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radiative collapse and developed dense post-shock layers (R), and shock that remained non- 
radiative (NR), are indicated. The right panels (Scenario 2) show the development of a shock 
(Vb = 3, 500 km s -1 , p HII = 150 amu cm -3 , d spot ~ 4 x 10 16 cm, p = 10 5 amu cm -3 ) where 
tcooi <C tcross- Due to the high density in this obstacle, all of the transmitted shocks undergo 
thermal collapse almost immediately 



-55 - 




E i i i ] 

45 90 135 180 

Obliquity [deg] 

Fig. 9. — The pressure behind the reflected shock as a function of obliquity immediately after 
impact of the blast wave for a spherical cloud having a density enhancement 5 = po/ pmi ~ 
70. We do not show the function for obliquities higher than 135° since the behavior at the 
back of a cloud is complex and highly time dependent due to the "healing" of the blast wave. 
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Fig. 10. — Range of transmitted shock velocities (V s h oc k) as a function of pre-shock density 
(p ) for three different estimates of the blast- wave pressure. The shaded grey region shows 
the range of transmitted shock velocities (due to obliquity) for our best estimate of the 
blast-wave pressure (V b = 3500 km s _1 , p H n = 150 amu cur 3 ). The upper boundary is the 
velocity of the transmitted shock at the tip of the protrusion and the lower boundary is the 
velocity at the side of the protrusion. The range of allowed shock velocities for our high 
blast-wave pressure estimate [Vb = 4100 km s _1 , pun — 250 amu cm -3 (dashed)], and the 
low estimate [Vb = 2800 km s" 1 , p HII = 100 amu cm -3 (dotted)} are also shown. 
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Fig. 11. — The top figure shows the temperature (thin) and density (thick) structure for 
electrons (dashed) and protons (solid) in a radiative shock of velocity V s = 250 km s _1 . The 
bottom figure shows the integrated line emissivities for N V A 1240 (dotted), C IV A 1550 
(long dashes), [N II] AA6548, 6584 (short dashes), and Ha (solid). The shading identifies the 
different regions of the shock as described in the text. 




Fig. 12. — Model shock cooling time t coo i versus shock velocity V s hock for initial obstacle 
density po = 2 x 10 4 amu cm~ 3 (solid). The power-law fit t coo i oc V^ 8 (dashed) is also shown. 
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Fig. 13. — Line production efficiencies rj for N V A1240 (dotted), C IV A1550 (long dashes), 
[N II] AA6548, 6584 (short dashes), and Ha (solid) as a function of shock velocity for fully 
cooled shocks (thick) and shocks which have only aged for 3 years (thin). 
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Fig. 14. — Similar to Figure 10, with additional lines (thick) dividing the phase space between 
radiative shocks and non-radiative shocks for cooling time t coo \ = 2 and 4 years. 
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Fig. 15.— Contours of equal line ratios for [0 I] (6300+6364)/5577 = 18.5 ±4.2 (long dash), 
[0 III] (4959+5007) /4363 = 6.6 ±3.4 (short dash), [N II] (6548±6584)/5755 = 2.0 ±0.2 (dot 
dash), [S II] (6717±6731)/(4068±4076) = 0.1 ±0.04 (dotted), and [S III] (9068±9530)/6312 
= 20.7 ± 10.9 (solid). The middle line corresponds to the determined ratio of the dereddened 
line fluxes while the other two curves bracket the la errors. Since the [O III] A4363 line may 
be blended with some Fe emission, the measured ratio is a lower limit. The arrow attached 
to the [O III] ratio contours shows how this may affect the derived temperature and density. 



-62 - 



NVA1240 



Si IVA1397 




2500 3000 3500 4000 4500 5000 

lime [days sines SN] 



ON] IU0D 




2500 3000 3500 4000 4500 5000 
lime [days sines SN] 

CIVJI1550 




2300 3000 3500 4000 4500 5000 
lime [days since SN] 




3000 3500 4000 4500 5000 

1ms [days sines SN] 



NIV] A1486 




3000 3500 4000 45 CO 5000 
time [days sines SN] 

He II A 1640 




3000 3500 4000 4500 5000 
time [days since SN] 



Fig. 16. — The observed (error bars) and modelled (Model 1: empty squares; Model 2: filled 
circles) light curves for the UV lines. The dashed line shows the earliest date (day 2932) 
that Spot 1 appeared as reported by Lawrence et al. (2000a). 
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Fig. 17. — Shock area growth for the V s = 135 km s -1 (filled sqaures) and V s = 250 km s _1 
(empty circles) shocks in the best-fit Model 2 scenario. The dashed line shows the earliest 
date (day 2932) that Spot 1 appeared as reported by Lawrence et al. (2000a). 
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Fig. 18. — Model shock-front geometries and flow fields for the radiative shocks in the 
obstacle (shock incident from the left side). The shock speed varies along the surface from 
250 km s _1 at the front to 100 km s~ x at the back. Cylindrical symmetry is assumed for all 
models. 
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Fig. 19. — Model Ha line profiles (solid) for the shock-front geometries in Figure 18 are 
shown. These profiles are smoothed with a Gaussian with FHWM = 100 km s _1 (dotted). 
The day 4568 Ha line profile from the medium resolution O'.'l G750M data is plotted for 
comparison (square). The uncertainties of the observed profile are smaller than the size of 
the symbol. The low- velocity region (—50 to +50 km s _1 ) is dominated by emission from 
the circumstellar ring (cf. Figure 4) and is therefore not displayed. 
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Fig. 20. — One-dimensional VH-1 simulation (800 spatial zones) of a radiative shock (V s ~ 
250 km s _1 ) developing in a dense (no = 3.3 x 10 4 cm' 3 ) obstacle. Shown is a space-time 
diagram where darker shade represents higher density. The overall flow is advecting to the 
left due to the inflow boundary conditions chosen for the right-side boundary (Vi n fi ow = 
220 km s" 1 ). 
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Fig. 21. — Two-dimensional VH-1 simulation (80 x 400 zones) of the development of a 
"turbulent" cold layer behind a radiative shock (V s ~ 250 km s _1 ) in a dense (no ~ 3.3 x 
10 4 cm -3 ) obstacle. The elapsed time between individual frames is one year. The obstacle 
has an imposed 5% sinusoidal density perturbation in a direction perpendicular to the shock 
front. 



